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ABSTRACT 

We study the small population of high-redshift (zgm > 2.7) quasars detected by GALEX, whose far-UV 
emission is not extinguished by intervening H I Lyman limit systems. These quasars are of particular impor- 
tance to detect intergalactic He II absorption along their sightlines. We correlate almost all verified Zem > 2.7 
quasars to the GALEX GR4 source catalog covering ~ 25000deg^, yielding 304 sources detected at S/N> 3. 
However, ^ 50% of these are only detected in the GALEX NUV band, signaling the truncation of the FUV 
flux by low-redshift optically thick Lyman limit systems. We exploit the GALEX UV color mpuv — 'Wnuv to 
cull the most promising targets for follow-up studies, with blue (red) GALEX colors indicating transparent 
(opaque) sightlines. Extensive Monte Carlo simulations indicate a He II detection rate of ^ 60% for quasars 
with mpuv — »Jnuv ^5 1 ^em ^ 3.5, a ~50% increase over GALEX searches that do not include color informa- 
tion. We regard 52 quasars detected at S/N> 3 to be most promising for HST follow-up, with an additional 1 14 
quasars if we consider S/N> 2 detections in the FUV. Combining the statistical properties of H I absorbers with 
the SDSS quasar luminosity function, we predict a large all-sky population of ^ 200 quasars with Zem > 2.7 
and / < 19 that should be detectable at the He II edge at m304 < 21. However, SDSS provides just half of the 
NUV-bright quasars that should have been detected by SDSS & GALEX. With mock quasar photometry we 
revise the SDSS quasar selection function, finding that SDSS systematically misses quasars with blue u—g<2 
colors at 3 < Zem ^3.5 due to overlap with the stellar locus in color space. Our color-dependent SDSS selection 
function naturally explains the inhomogeneous u — g color distribution of SDSS DR7 quasars as a function of 
redshift and the color difference between color-selected and radio-selected SDSS quasars. Moreover, it yields 
excellent agreement between the observed and the predicted number of GALEX UV-bright SDSS quasars. We 
confirm our previous claims that SDSS preferentially selects 3 < Zem ^3.5 quasars with intervening H I Lyman 
limit systems. Our results imply that broadband optical color surveys for 3 < Zem ^3.5 quasars have likely 
underestimated their space density by selecting IGM sightlines with an excess of strong H I absorbers. 
Subject headings: diffuse radiation — intergalactic medium — quasars: absorption fines — surveys — tech- 
niques: photometric — ultraviolet: galaxies 



1. INTRODUCTION 

The intergalactic space is pervaded by a filamentary cos- 
mic web of gas of almost primordial composition, the so- 
called intergalactic m edium (IGM) , seen in absorpt ion against 
background sources (lRauchlll998l: lMeiksinll2009l) . The ab- 
sence of H I Lya absorption troughs in spectra of Zem < 6 
quasars signals that the h ydrogen in the IGM is highly ionized 
(iGunn & Petersonl[i965i) . Instead, the plethora of narrow H I 
Lya absorption lines, known as the Lya forest, traces the tiny 
residual neutral hydrogen fraction of the IGM as the largest 
reservoir of baryons in the universe. The ionizing radiation of 
quasars and star-forming galaxies is filtered by the IGM, lead- 
ing to the buildup of the UV background radiation field that 
determines the ionization state of the gas (Haardt & Madau 
[l996; Fai-daletal. 1998; Faucher-Giguere et alJUOOSjT The 
UV background changes in amplitude and spectral shape due 
to evolution in the source number d ensity, cosrnologic al ex- 
pansion and structure formation (e.g. 'D ave et al.lll999b . This 
is particularly important for the ionization state of helium, 
the second most abundant element in the IGM. Due to its 5.4 
times higher recombination rate and 4 times higher ionization 
threshold, the reionization epoch of helium (He II — >He III) 
is expected to be delayed with respect to hydrogen. 

The Lya transition of intergalactic He II at Aiest = 303. 78A 
is observable in the far UV (FUV) from space only at z > 2 
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due to the Galactic Lyman limit. The determination of the 
He II reionization epoch via the He II Gunn-Peterson test 
towards high-redshift quasars has been a major goal in ex- 
tragalactic UV astronomy since the laun ch of the Hubble 
Space Telescope (HST , e.p. lMiralda-Escude & Ostr iker 199(J; 
iMiralda-Escudel [l993h . However, the accumulated Lyman 
continuum (LyC) absorption of the H I absorber population 
severely attenuates the quasar flux in the FUV, rendering just 
a fe w percent of Zem > 3 sightlines to be relatively transpar- 
ent ('M 0ller & Jak obsen* 1990). The combination of the rising 
LyC absorption and the declining quasar luminosity function 
results in a sharply dropping number of observable UV -bright 
quasars at Zem > 3 jPicard & Jakobsenlll993l:ljakobseniri998,) . 

Until very recently He II Lya absorption had been found 
only in a handful of sightlines despite considerable effort, 
since the UV fluxes of most targeted quasars had been un- 
known. HST observations of Q 0302-003 at Zem = 3.285 
(iJakobsen et aH [19941: iHoean et alJ [19971: iHeap et alJ llOOOh 
and PKS 1935-692 at z.n. = 3.18 (lAnderson et al.lll999h re- 
vealed a high He II effective optical depth at z > 3 that is con- 
sistent with a Gunn-Peterson trough (Teff.Heii > 3). In con- 
trast, the lines of sight towards HS 1 700+6416 at Zem = 2.736 
(iDavidsen et lI\ \T99^ iFechner et aL 2006), HE 2347 -4342 



at Zem = 2 885 dReimers et alj 119971: iKrisset ai] I2OOII: 
iSmette et al.ll2002l: iZheng et alJl2004bl: iShuU et alJl 2004') and 
HS 1157+3143 at Zem = 2.989 (iReimers et all 12005) show 
patchy He II absorption with voids (leSMeii < 1) ™d troughs 
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(Teff Heii > 3). At z < 2.7 this patchy absorption evolves 
into a He II Lya forest that has been resolved in high- 
resolution spectra obtained with the Far Ultraviolet Spectro- 
scopic Explorer ( FUSE, Kriss et al. 2001; Zheng et al. 2004^ 
IShull et al.ll20(Mt iFechner et al. 2006). 

The strong evolution of the He II absorption suggests a 
late reionization epoch of helium at z ~ 3, when quasars 
have been sufficiently abundant to supply the required hard 
photons. The patch-work of absorption and transmission 
evokes a picture of overlapping He III zones around quasars 
that lie close to the sightline (Reimers et al. 1997; Heap et al] 
120001; ISmette et all l2002h . Indeed, the He III proximity 
zones of quasars hav e been detected both along the line 
of sight dHogan et alj [19 97; Anderso n et all 119991) and in 
transverse direction Jjakobsen et ahi2003h . In the past few 
years, great progress has been made in developing the the- 
oretical framewor k to interpre t these observations. Both 
semi-analytic (e. g.l Haardt & M adau 1996; Fardal et al. 1998; 
iGleser et al.l 12005; Furlanetto & Dixon 2010) and numeri- 
cal radiative tra nsfer si mulations ( | Ma selli & Ferrara 2005; 
Tittlev & Meiks in 2007; Paschos et alj |2007; McOuinn et al. 
2001 indicate that the He II reionization process should be 
very inhomogeneous and extended over 3 < z < 4, since 
rare luminous quasars dominate the photoionizing budget of 
the overall quasar population. The few quasars contribut- 
ing to the UV radiation field at the He II ionization edge 
at a given point likely give rise to fluctuations in the FUV 
background that can be tracked by the co-spatial absorption 
of He II a nd H I dB oltonetaT '2006'; 'Wor seck & Wisotzkil 
12006; Worseck et a ll ,2007 ; Furlanetto 2009). Th e UV back- 
ground hardens as He II reionization proceeds dHeap et alj 
l2000t IZheng etani2004bl) . but > 10 Mpc fluctuations are ex- 
pected to persist even after its end (IFechner & Reimersl2007h . 

Other, more indirect observations might suggest that He II 
reionization is ending at z ^ 3. The IGM is reheated as the 
individual He III bubbles around quasars overlap, however 
the amplitude of thi s t emperature jurn p is highly uncertain 
dBolton et alj l2009allbl ; iMcOuinn et al] [2009). Observation- 
aUy, several studies indicated a jump in the IGM temperature 
atz ~ 3 (Hicotti et al. 2000; Schaye et al. 2000; Theuns et al., 
I2002h . whereas others are consistent with an almos t con- 
stant IGM temperature at 2 < z < 4 (M cDonald et al.l 1200 ll: 
iLidz et an r2010). Moreover, photoionization models of metal 
line systems indic ate a significant hardening of t he UV back- 
ground at z < 3 dAgafonova et all I2005L l2007h . However, 
these observations are restricted to rare metal line systems 
showing various ions with a simple velocity structure. 

At present, the five He II absorption sightlines studied at 
scientifically useful spectral resolution provide the best ob- 
servational constraints on He II reionization. However, just 
one or two sightlines probe the same redshift range, and given 
the large predicted variance in the He II absorption, this small 
sample clearly limits our current understanding of He II reion- 
izatioifl. The Sloan Digital Sky Survey (SDSS) has dra- 
matically increased the number of high-redshift quasars to 
search for the presence of flux at He II Lya, yielding three 
Zem > 3.5 quasars w ith detected He II Gunn-Peterson troughs 
dZheng et al."2004a', lOOS, 2008). More importantly, the al- 
most completed first UV all-sky survey with the Galaxy Evo- 
lution Explorer (GALEX) enables the pre-selection of UV- 
bright quasars for follow-up UV spectroscopy, leading to the 
recent discovery of 22 new clear sightlines towards SDSS 

' Ironically, the z ~ 6 epoch has substantially better statistics. 



quasars at 3.1 < Zem < 3.9 dSvphers et alj2009allH) . The avail- 
able GALEX photometry dramatically increases the survey 
efficiency by almost an order of magnitude to ~ 42% in the 
ISyphers et al.l survey. 

The recently installed Cosmic Origins Spectrograph (COS) 
on HST offers unprecedented sensitivity to study He II reion- 
ization via He II Lya absorption spectra. With its confirmed 
throughput at A > 1 105 A ( iMcCandhss etalj2010l) HST/COS 
is now able to probe He II Lya at z > 2.64, thereby covering 
the full redshift range of int erest for He II reionization. Very 
recently, [ShuU et al, d2010l) presented a high-quality COS 
spectrum of HE 2347—4342, dramatically improving on ear- 
lier FUSE data. In the near future, COS will be employed 
to both obtain follow-up spectroscopy of the recently con- 
firmed He II sightlines, and to discover new ones. In this paper 
we introduce the quasar UV color measured by GALEX as 
a powerful discriminator to select the most promising sight- 
lines for follow-up spectroscopy. Moreover, we significantly 
improve on earlier predictions on the number of UV-bright 
quasars dPicard & Jakobsenll 19931: iJakobsenll 19981) . based on 
observational advances to characterize both the quasar lumi- 
nosity function and the optically thick IGM absorber distribu- 
tion. The structure of the paper is as follows: In [j2]we will 
present our sample of verified high-redshift quasars detected 
by GALEX. Section [3] describes our Monte Carlo routine to 
compute H I absorption spectra and to perform mock GALEX 
and SDSS photometry. In |j4]we determine the expected num- 
ber of UV-bright Zem > 2.7 quasars and establish GALEX UV 
color selection criteria to select quasars with probable He II- 
transparent sightlines. We compare the observed and pre- 
dicted number counts of UV-bright SDSS quasars in i|5]before 
concluding in ^ 

2. OUR SAMPLE OF Z > 2.7 QUASARS DETECTED BY 

GALEX 

2.1. The initial quasar sample 

We compiled a list of practicaUy all known quasars at 
Zem > 2.7 from four quasar samples . We s tarted with 
the SDSS DR5 quasar catalog (Schnei der et all 2007) and 
added all o ther spectroscopic SDSS targets from DR6 
dAdelman-McCarthv et al. 2008) and DR7 ( Abazaiian et al] 
12001 " identified as Zem > 2.7 quasars by the SDSS spec- 
tro Id pipeline. We supple mented this SDSS q uasar list by all 
Zem > 2.7 sources from the lVeron-Cettv & Ve ron (2006) cat- 
alog not discovered or verified by SDSS. This merged quasar 
catalog is inhomogeneous due to several reasons: (i) the 
SDSS DR5 quasar catalog represents a non-statistical sample 
due to changes in the quasar selection criteria in the course 
of the SDSS dRichards et alJ 120061: iSchneider et all [2007), 
(ii) the inclusion of SDSS quasars discovered by serendip- 
ity (Stoughtonet al. 2002), (iii) the redshifts of most SDSS 
DR6/7 sources have not been verified by eye, and (iv) the 
iVeron-Cettv & Veron (2006) catalog is inherently inhomoge- 
neous as it is a collection of quasars discovered by various 
surveys with sometimes unknown selection criteria. 

The merged list of quasars contained 12373 unique en- 
tries. However, among them there are SDSS DR6/7 sources 
misidentified as high-z quasars by the SDSS source identifica- 
tion algorithm either due to misclassification or a wrong red- 
shift assignment. We refrained from the tedi ous visual clas- 
sification of all spectrold DR6/7 quasars (see ISchneider et al] 
120 lOl for the DR7 quasar catalog compiled after our analysis 
was finished), and limited our visual verification to the sub- 
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set of SDSS DR6/7 sources actually de tected by GALEX (see 
below). Moreover, we caution that the I Veron-Cettv & VeronI 
( 120061) catalog contains a fair number of quasar candidates 
with estimated redshifts from slitless spectroscopic surveys. 
Many of these redshifts will be grossly overestimated as most 
slitless spectroscopic surveys assign the highest plausible red- 
shifts if just a single emission line is present. Consequently, 
we removed all misidentified SDSS sources and all quasar 
candidates without unambiguous redshifts from follow-up 
spectroscopy, but only after cross-correlating the initial quasar 
sample to the GALEX GR4 source catalog. 

2.2. Cross-correlation with GALEX GR4 

The GALEX satellite currently performs the first large- 
scale UV imaging survey (Ma rtin et al.ll2005l iMorrissey et alJ 
|2007). Most images are taken simultaneously in two broad 
bands, the near UV (NUV, -1770-2830A) and the far UV 
(FUV, -1350-1780A) at a resolution of - 5" full width at 
half maximum (FWHM). Three nested GALEX imaging sur- 
veys have been defined: the All-Sky Survey (AIS) cover- 
ing essentially the whole extragalactic sky (^^26000 deg^) 
to mAB '^21, the Medium Imaging Survey (MIS) reaching 
'wab ^ 23 on 1000 deg^, and the Deep Imaging Survey (DIS) 
extending to hjab ^ 25 on 80 deg^. These main surveys are 
complemented by guest investigator programs. The GALEX 
Data Release 4 (GR4) covers -25000 deg^, 96% of the antic- 
ipated AIS survey area. The officially distributed GR4 data 
has been homogeneously reduced and analyzed by a dedi- 
cated software pipeline. A previous version of this pipeline 
used for the earlier GR 3 data release is described in detail by 
IMorrissey et alJ (l2007h . 

We cross-correlated our initial quasar list to the available 
GALEX GR4 source catalogs using a maximum match radius 
of 4.8" around the optical quasar position. The match radius 
approximately corresponds to the typical GALEX FWHM 
and was chosen to account for the degrading astrometric ac- 
curacy of GALEX towards the detection limit where we ex- 
pect most of the rare UV-transparent quasars (see fj23] be- 
low). In comparison, the pos itional error s of the quasars are 
negligible, 0"! for SDSS (iPier et al.ll2003l) and < 1" for the 
|Veron-Cettv & Veron catalog quasars. 

2.3. Source verification and catalog completeness 

Substantial screening of the cross-matches was required 
to create our final list of real Zem > 2.7 quasars detected in 
GALEX GR4. We visually confirmed the redshift of ev- 
ery detected SDSS so urce and searched the references of the 
IVeron-Cettv & V6roni catalog quasars for unambiguous red- 
shift determinations a nd plotted spectra. A large fraction of 
the GALEX-detected IVeron-Cettv & VeronI quasars had un- 
confirmed slitless spectroscopic redshifts, in line with our as- 
sertion that most of them are in fact low-redshift interlop- 
ers. Consequently we removed these unconfirmed candidates. 
In addition, we flagged obvious broad-absorption-line (BAL) 
quasars which are rarely usable for IGM studies due to the dif- 
ficulty in disentangling the IGM absorption along their sight- 
lines from the high-velocity quasar outflows. This flagging 
was somewhat restrictive, as it was based on the visual ap- 
pearance of the spectrum (if available), and quasars with con- 
fined low-velocity narrow BAL systems were kept in the sam- 
ple. FinaUy, we inspected the SDSS images of all GALEX- 
detected quasars in the SDSS DR7 footprint, and flagged 
cases of potential source confusion with blue optical neigh- 



bors at < 5" separation caused by the broad GALEX point 
spread function (PSF). Specifically, a quasar was flagged if 
the spectral energy distribution of the neighbor (as estimated 
from the SDSS photometry) was likely to extend to the UV 
(e.g. significant u band flux). In total, ~20% of the SDSS 
quasars were flagged. La cking deep multi-band photometry, 
we could not inspect the IVeron-Cettv & VeronI quasars out- 
side of the SDSS footprint with the same scrutiny. For quasars 
imaged in multiple GALEX exposures we kept only the most 
significant detection, usually in the deepest exposure unless 
affected by obvious image artifacts. For every source for- 
mally detected in only one GALEX band we obtained a 1 C7 
upper limit on the flux in the other In total, we were left 
with 803 verified Zem > 2.7 quasars with Ukely GALEX GR4 
counterparts. Almost all of them (782) have been imaged in 
both GALEX filters, allowing for constraints on the UV color 

Due to the strong Lyman continuum absorption by the in- 
tervening IGM most of these high-redshift quas ars a re faint 
in the UV even if they are optically bright (see i l4. II below). 
Most of these rare high-redshift quasars with appreciable UV 
flux will be detected at low signal-to-noise (S/N) close to 
the limits of the defined GALEX imaging surveys. Incom- 
pleteness arises in the source catalog at low S/N, resulting 
in false negatives (nondetections in one or both bands) and 
false positives (no UV flux at all). The low-S/N UV fluxes 
are naturally uncertain and li kely overestimated due to Ed- 
dington bias (IMorrissey et al.l l2b07). The detection repeata- 
bility is generally low at the survey limit, and the detectabil- 
ity of sources sometimes depends on subtle c hanges in the 
data analysis. For example, two quasars that ISvphers et alj 
(l2()09b.) confirmed to show flux at He II Lya were listed in 
the GRI catalog, but not in further GALEX data releases with 
improvements in survey depth, calibration and source detec- 
tion routines. While low-S/N detections might still indicate 
UV-transparent quasars, we limit our statistical studies (!|5]l to 
sources with S/N> 5 in at least one of the GALEX bands. At 
the lowest S/N ratios encountered one has to question the real- 
ity of the UV detection, in particular if a source is seen just in 
one GALEX band. Sources formally detected in both bands 
should be less affected, as source detection is performed in- 
depen dently on the FUV and NUV images (IMorrissey et al.l 
l2007h . Compared to the general incompleteness at faint mag- 
nitudes, the subtle effect of PSF and sensitivity degradation 
at the rim of the GALEX field of view can be neglected. 
We therefore performed our correlation analysis on the full 
GALEX tiles, thereby maximizing the number of promising 
UV-bright quasars for He II studies. 

We investigated the astrometric performance of GALEX 
in the low S/N regime by calculating the offset between the 
optical quasar catalog position and the GALEX NUV and/or 
FUV position. Given the nested GALEX surveys with a large 
spread in depth, the astrometric accuracy primarily depends 
on S/N rather than on magnitude. Figure [T] plots the cumu- 
lative fraction of the squared separation between the GALEX 
positions and the optical position of GALEX-detected SDSS 
Zem > 2.7 quasars for various ranges in S/N. In this metric, 
false positives will be uniformly distributed in r^, whereas 
quasar (neighbor) matches should be concentrated at small 
(large) offsets. Indeed, for SDSS quasars having blue opti- 
cal neighbors within 5", the distribution has two peaks, one 
at small separations for matches to the quasar, and one at 
large separations corresponding to the detected blue neigh- 
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Figure 1. Cumulative fraction of the squared separation between the 
GALEX positions and the optical position of GALEX-detected SDSS 
quasars. Thick (thin) black lines show FUV-optical (NUV-optical) distribu- 
tions for various ranges in S/N. The thin gray line shows the cumulative distri- 
bution of the NUV-optical separations of SDSS quasars having blue optical 
neighbors within ~5". The diagonal line denotes the uniform distribution 
with squared separation that is expected for false positives. 



bor instead of the quasar. Therefore, it is essential to flag 
such cases of potential source confusion caused by the broad 
GALEX PSF. With the assumption that all GALEX sources 
in the SDSS footprint should have SDSS counterparts, the 
GALEX sources without sufficiently blue optical neighbors 
are either UV counterparts to the quasars in our catalog or 
false positives (noise). 

Figure [T] shows that for SDSS quasars without blue optical 
neighbors the distributions peak at small offsets with a clear 
dependence on S/N. Almost all FUV (NUV) S/N> 5 detec- 
tions are within r < 2" (r < 3") of the optical position with 
the difference being du e to the better resolution in the FUV 
jMorrissev et al.ll2007h . At lower S/N the astrometric accu- 
racy degrades and the rate of false positives should increase. 
At S/Nnuv < 2 the cumulative fraction begins to resemble 
the one expected for false positives, with the excess indicat- 
ing some real detections among them. Since the offset distri- 
butions at S/Nnuv > 2 are much more concentrated, we infer 
that a limiting S/N> 2 rather than a fixed limit in the matching 
radius yields a source catalog of high purity and complete- 
ness. Our chosen matching radius of 4.8" likely encompasses 
all true matches with S/N> 3, whereas a few real 2 <S/N< 3 
detections (without neighbors) might exist at even larger sep- 
arations. After excluding 117 (~ 20%) of the SDSS quasars 
with neighbors, restricting our catalog to S/N> 2 (S/N> 3) 
in at least one GALEX band reduces the number of potential 
(probable) detections to 601 (304). 

We examined the GALEX source counts within 3' around 
our quasars to estimate the probability of residual false 
matches between quasars and GALEX detections. Despite 
their low resolution, GALEX images are confusion-lim ited 
only in the longest DIS exposures (Ha mmer et al.ll20Tol) due 
to the low source density in the UV. The measured den- 
sity of S/Npuv > 2 detections in a typical MIS exposure 
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Figure 2. Lick/Kast spectrum of the emission line galaxy J 1943— 1502 (z = 
0.0427). Identified emission lines are marked. 

is ^ 1/arcmin^, which accounts for both real sources Q and 
false positives. At this low of a source density, the chance 
for any S/Npuv > 2 detection to fall in our 4.8" aperture is 
small (< 2%). Given that the source density on AIS plates is 
even lower, we conclude that essentially all FUV matches on 
AIS and MIS plates will correspond to optical sources within 
the chosen aperture. The rejection of SDSS quasars with 
blue neighbors probably excluded several real SDSS quasar 
matches (Fig.[TJ, so that we consider > 98% of the remaining 
FUV-SDSS matches to be real. For non-SDSS quasars the re- 
maining source confusion is more important than the rate of 
spurious detections. Adopting our SDSS neighbor fraction of 
^ 20%, we estimate a purity of ~ 80% for the quasars not 
imaged by SDSS. Due to the challenging reduction and anal- 
ysis of DIS plates, we flagged the 23 quasars detected on DIS 
plates as still potentially affected by source confusion (only 7 
are in the constrained sample discussed in ^A.2\ . 



2.4. Comparison to Source Matching in \Syphers et al\ 
i.2009a) 

Recentlv. lSvphers et al.l (l2009ah pubhshed a catalog of 593 
sources detected in GALEX GR4 and its small extension 
GR5. Apart from a slightly higher redshift cutoff (z > 2.78) 
and a smaller matching radius (3" around the quasar), their ap- 
proach to source matching (not target selection) was similar to 
ours. However, they admitted that they did not verify the red- 
shifts of the 165 sources with GALEX GR4-I-5 counterparts 
stemm ing fro m the Veron-Cettv & Veron (2006) catalog. 
ISvphers et all f2009al) presented follow-up HST/ACS UV 
prism spectroscopy of one of these, J 1943 — 1 502, with an es- 
timated slitless spectroscopic redshift of 3.3 dCrampton et al.l 
|1997|) . In order to establish whether this object can be used 
for He II IGM studies, we obtained an optical spectrum with 
the Kast spectrograph at the 3-m Shane Telescope at Lick 
Observatory. We confirm J 1943 — 1502 as a naturally UV- 
bright low-redshift emission li ne galaxy rather than a quasar 
(Fig, in. We caution that th e ISvphers et all ( l2009ah Ust of 
IVeron-Cettv & VeronI ( 120061) sources contains 41 more such 
candidates the redshifts of which should be confirmed before 
embarking on follow-up UV spectroscopy with HST. In ad- 
dition, 5 other sources from the Veron-Cettv & Veron (200^ 
catalog that are listed by Svphers et al.. (2009a) as GALEX- 
detected Zem > 2.7 quasars are actually at lower redshifts ac- 
cording to our visual inspection of their spectra. 

^ We compared our measured source density to the literature (Bianchi et ajj 
1^07; Hammer et al. 2010). At our low S/N threshold we only recover ~ 60% 
of the predicted sources on a given GALEX plate due to incompleteness at 
the survey limit. 
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Figure 3. HST spectra of 12 quasars (gray) and NUV & FUV fluxes measured by GALEX (filled circles with eiTor bars and indicated bandpass). Aitows indicate 
upper limits from GALEX non-detections. Dashed lines mark zero flux. The thick vertical lines indicate the expected onset of He II absorption. 



2.5. The GALEX UV colors of high-redshift quasars 

The large sky coverage of GALEX enables the recovery of 
many UV-bright Zem > 2.7 quasars that have previously been 
followed up with HST to search for He II absorption by the 
IGM. GALEX recovers all 8 quasars known to show flux 
at He II Lya that had been selected for observations before 
the launch of GALEX. S vphers et al.l (|2009a b) recently con- 
firmed 22 GALEX- selected sighthnes to show He II, and all 
but the two listed only in GRl are contained in the GALEX 
GR4 source catalog. 1 3 of the total 30 confirmed He II quasars 
are detected by GALEX at a low S/N< 3, and we suspect that 
there is a larger population of UV-transparent quasars missed 
at the GALEX survey limit. We also recovered UV-bright 
quasars considered in previous photometric and spectroscopic 
surveys for He II with HST, the sightlines of which are inter- 
cepted by optically thick Lyman limit systems redward of the 
onset of He II absorption. 

In Fig. [3] we compare the GALEX fluxes of 12 quasars to 



their UV spectra taken with HST. Their GALEX UV magni- 
tudes are provided in Table [T] together with references to the 
UV spectra and the Lyman limit systems zeroing the spectral 
flux (if any). As these quasars are bright in the UV they are 
imaged with GALEX at high S/N, so that the GALEX fluxes 
are in very good agreement with the HST spectrophotome- 
try. More interestingly, we find that several opaque sight- 
lines are just detected in the NUV, but not in the FUV as 
expected (LBQS 0041-2603, Q 0055-269 and UM 366 in 
Fig EJl. In contrast, quasars that show flux down to the onset 
of He II absorption are detected in both bands with the flux 
rising towards shorter wavelengths as it recovers from par- 
tial Lyman Hmit systems (HS 1700+6416 and Q 0302-003). 
Thus, the GALEX UV color mpuv — wnuv can be used to effi- 
ciently distinguish between opaque sightlines (red UV color) 
and transparent ones (blue UV color). The only quasars 
that remain insensitive to this obvious color selection crite- 
rion are those caught by an optically thick Lyman limit break 
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Table 1 

Data on the UV-bright quasars shown in Fig. [3] 



Object 




mpuv [AB] 


'"Nuv [AB] 


HST spectram 


ZLLS 


References 


PKS 2212-299 


2.706 


20.74 


20.35 


STIS G230L 


0.6329 


Rao et al. (2006) 


HS 1700+6416 


2.736 


18.94 


18.74 


FOS G130H/G190H/G270H 




Reimers et al. (1992); Evans & Koratkar (2004) 


UM682 


2.756 


19.68 


19.65 


FOS G160L/PRISM 




HST Archive 


Q 0903+175 


2.773 


> 23.32 


21.50 


FOS G160L/G270H 


BAL 


Tumshek et al. (1996_) 


Q 0207-398 


2.805 


20.29 


21.07 


FOS G160L/G270H 




Bechtold et al. (2002) 


LBQS 0041-2638 


3.053 


> 24.32 


21.67 


STIS G230L 


1.38 : 


HST Archive 


UM366 


3.141 


> 24.40 


21.35 


FOS G160L/G270H 


1.6128 


Rao & Tumshek (2000): Evans & Koratkar (2004) 


HS 1140+3508 


3.147 


20.68 


19.92 


STIS G140L 


0.557 


HST Archive 


UM670 


3.163 


20.97 


20.22 


FOS G160L 


0.47 : 


Lvons et al. (1994): Evans & Koratkar (2004) 


Q 0302-003 


3.285 


21.37 


21.76 


STIS G140L/G230L 




Jakobsen et al. (1994); Heap et al. (2000) 


PKS 1442+101 


3.530 


20.98 


20.26 


FOS G160L/PRISM 


0.621 : 


Lvons et al. (1995); Evans & Koratkar (2004) 


Q 0055-269 


3.665 


> 23.17 


20.97 


FOS G160L/PRISM 


1.5335 


Cristiani et al. (1995); Evans & Koratkar (2004) 



just in the narrow range between the GALEX FUV band 
and the onset of He II absorption (HS 1 140+3508, UM 670, 
PKS 1442+101, PKS 2212-299 in Fig. |3]i. We also iden- 
tify two FUV-detected quasars, the HST spectra of which do 
not extend to He II Lya in the rest frame of the quasar, lo- 
cated near the UV sensitivity cutoff of HST (UM 682 and 
Q 0207—398). These two sightlines are likely transparent, as 
there are no obvious strong Lya absorbers that could cause 
a Lyman limit break in the 200A gap to the onset of He II 
absorption. 

With the additional quasars targeted in recent surveys for 
He II sightlines (Syphers et al. 2009a b) we can confirm the 
trend that most quasars with flux down to He II Lya show blue 
GALEX colors, whereas most fruitlessly targeted quasars are 
characterized by red colors (see Fig. [12] below). Although 
more uncertain at low S/N, the colors still distinguish both 
quasar populations at S/N> 3. Excluding sources with neigh- 
bors, ^ 50% of the SDSS quasars in our sample are detected 
at S/N> 3 in the NUV band, but are lacking a significant FUV 
detection (S/Npuv < 2), indicating the ubiquitous strong Ly- 
man continuum absorption. In particular, FUV dropouts de- 
tected in the NUV at high significance likely correspond to 
optically thick Lyman limit breaks. 

In the following sections we will further explore how to fur- 
ther constrain our sample by the GALEX UV color to select 
the most promising quasar sightlines to detect He II absorp- 
tion. This requires one to create mock quasar spectra with 
appropriate H I absorption, and to perform GALEX photom- 
etry on them to relate the GALEX UV color to the Lyman 
continuum absorption along the line of sight. 

3. MONTE CARLO SIMULATIONS OF 
HIGH-REDSHIFT QUASAR SPECTRA 

3.1. Monte Carlo model for the H I Lyman series and 
Lyman continuum absorption 

3.1.1. General procedure 

For the problem at hand we followed standard practice to 
generate Monte Carlo (MC) H I Lyman forest and Lyman con- 
tinuum absorption spectra from the observed statistical prop- 
erties of the Lya forest (e.g . M0ller & Jakobsen 1990; Madau 
[1991 iBershadv et an[T999i: ilnoue & Iwata 2008). The spec- 
tra were generated under the null hypothesis that the Lya for- 
est can be approximated as a random collection of absorption 
lines (Voigt profiles) with uncorrelated parameters (redshift z, 
column density A^hi and Doppler parameter b). From the line 
list representing the H I absorber population on a given line 
of sight from z = to an emission redshift Zem we created ab- 



sorption spectra of the Lyman series (up to Ly30). Individual 
resolved Voigt profiles were computed on AA = O.OSA pix- 
els using the approximation by Tepper-Garcia ( 2006). Lyman 
continuum abso rption was included using the H I ionization 
cross section by Verner et a l. ( 1996). 

In order to accurately predict the far-UV attenuation of 
high-redshift quasars by the IGM we desired a model that 
successfully reproduces the observed statistical properties of 
the Lya forest at all redshifts, in particular concerning high- 
column density absorbers. Considering the recent observa- 
tional advances in Lya forest statistics, we deviated from pre- 
vious simple MC descriptions of the Lya forest and adjusted 
our input parameters as detailed in the following. 

3. 1 .2. The absorber redshift distribution function 

In our MC model the number of H I absorbers per line 
of sight in a given redshift range is a Poisson process 
(Zuo & Phinnev 1993). The observed mean differential line 
density per unit redshift is commonly parameterized as a 
power law dn/dz| forest °^ (1 +z)^ that results in an effective 
optical depth Tf^ g °^ (1 + zF ^' for Lya (and higher order 
series) absorption dZuoll 19931) . While there is some evidence 
that the redshift evolution depends on the column density even 
in the low-column density Lya forest, the uncertainties are 
still large due to the non-unique process to deblend the forest 
into a series of Voigt profiles especially at z > 3, incomplete- 
ness at the lowest column densities (logA^Hi 12.5), and the 
paucity of moder ate-co lumn density (logA^ni ^ 14.5) systems 
(IKim et al.iil99'7l [20021 ). We therefore chose to parameterize 
dn/dz for absorbers with 12 < logA^ni < 19 as a single power 
law, the parameters of which were fixed by requiring each 
simulated spectrum to be consistent with a specified power 
law in Teff_a(z). Observations point to a break at z ^ 1.5, be- 
low which there is litt l e evolution both in the line density (e.g . 
IWevmann et ai][T998t iKim et alj[2002l: iJanknecht et ani200a) 
and the mean a bsorption in the Lya forest Da = 1 — e^^'*" 
(iKirkman et al.l [2007). Thus, we assumed a broken power 
law for Teff,a(z). Knowing that a power-law line distribu- 
tio n generally will not y ield a power law for £>a(z) assumed 
by IKirkman et al.1 (EMzI), we converted their Da to Tgff a 

and obtained a fit Teffa(z) = 0.017(1 +z)'-^'' for z < 1.6. 
At 2 < z < 4 Teff,o; has_been preci s ely measured in high- 
resolu tion spectra dK im et al.' I2007t iFaucher-Giguere et al.l 
.2008; iDaU'Aglio et al. 2008), and the remaining disagree- 
ment at z > 4 is likely due to continuum uncertainties, where 
very few pixels remain unabsorbed even in high resolution 
spectra. We adopted the fit Teff,c( = 0.0062(1 +z)^ '''^ from 



GALEX far-UV color selection of UV-bright high-redshift quasars 



7 



iDair AgUo etal] ( 120081) . valid at 1 .8 < z < 4.6. Note that the 
break redshift cannot be determined as the intersection of the 
two power laws, since this would require one to extrapolate 
feff.aCz) beyond the quoted validity ranges. Since the break 
is observationally not well constrained, given the large scatter 
of Teff a measur ements at 1 .7 < z < 2 and the paucity of data 
at z 1.5 (Kir kman et alJl2007l) . we adopted a break redshift 
of z = 1 .5 for the broken power law in Teff a (z)- 

For logA^Hi > 19 absorbers we had to assume different red- 
shift evolution laws, both because these systems are gen- 
erally excluded in fits of Teff,a(z), and due to the fact that 
their number densities seem to evolve much slower with red- 
shift. For damped Lya systems (DLAs, logA^Hi > 20.3) 
we adopted dw/ dzlni.A = 0.044(1 +z)''^^, determined by 



iRao et alJ (l200g) over the redshift range < z < 5. Figure |4] 
co mpares the observed number densities of DLAs compiled 
by Ra o et all (|2006) to mock number densities obtained on 
4000 MC sightlines assuming their fit for dn/dzloLA- For 
Super Lyman Limit systems (SLLSs, 19 < logA^ni < 20.3) 
there are significantly less constraints in the literature. A 
maximum-likelihood power-law fit to the SLLS survey by 

lO'Meara etal] (12007) yields dn/dz|sLLS = 0.034(1 +zf-^'^ at 
1.8 < z < 4.2, but extrapolation to lower redshifts underesti- 
mates the lower limit nsLLS ^ 2nDLA at z < 1 .65 given by 
iRao et all (12006 ). Rather than a break in the number den- 
sity of SLLSs, this probably indicates that a much larger red- 
shift range is needed to accurately describe the number den- 
sity evolution of the rare SLLSs. By constraining the slope 
to 1.27 < 7sLLS < 2.14 (i.e. between the evolution rate of 
DLAs and the SLLS fit at high z), and considering the esti- 
mated total number of z < 1.65 SLLSs by Rao et al. (200^ 
we obtained a rough constraint on the low-redshift evolu- 
tion of SLLSs (dotted lines in Fig. [4 l ). Af ter binning the 
high-z measurements bv lO'Meara et air(l2007h we determined 

dn/dz|sLLS - 0.066(1 +z)' ™ by eye (Fig. gli, noting that 
these numbers are quite uncertain as the z < 2 SLLS popu- 
lation is not well constrained. 

3.1.3. The Doppler parameter distribution function 

Although the Doppler parameter distribution function 
d« / db is not required to calculate the attenuation of quasars 
by the IGM below the Lyman limit, our MC simulations re- 
produce the observed effective optical depth in the Lya for- 
est instead of a line density distribution. As the equivalent 
width of the lines on the flat part of the curve of growth 
(13.5 < logA^Hi ^ 18.5) depends both on the column den- 
sity A'hi and the Doppler parameter b, the line density that 
is consistent with our adopted Teff c((z) implicitly depends 
on the Doppler parameter distribution. For simplicity, we 
adopted the single parameter distribution function suggested 
by Hui & Rutledge (1999), dn/db oc b-^exp {-b'^/b^), with 
ba — 24kms^' dKim et alj 120011) independent of redshift 
and column density, and restricted to the plausible range 
lOkms'i <b< lOOkms-i. 

3. 1 .4. The column density distribution function 

Previous studies on the IGM attenuation of high-redshift 
sources approximated the column density distribution func- 
tion (CDDF) by a single or a broken power law, mainly 
driven by the reasonable approximation of the CDDF as 

a single power law dn/d^Vni °^ N^i with )3 ~ 1. 5 over 
practically the full observable column density range dTvtled 



0.1 



• SLLSs (obs.) 
■O SLLSs 
■ DLAs 
□ DLAs 




Figure 4. Adopted differential number densities dn/dz of SLLSs and DLAs 
as a function of redshift z- Filled symbols represent observed data (see text), 
whereas open symbols show the distributions recovered from 4000 MC sim- 
ulations. The straight solid and dashed lines denote the power law fits to the 
observed data adopted for the simulations. The dotted lines show different 
normalizations to yield the number density of SLLSs at z < 1 .65 observed by 
Rao et al. ( 2006) adopting the slope of the DLA evolution law and the SLLS 
slope found at high z. 



119871) . However, more recent studies revealed significant 
deviations in the high-redshift CDDF from a single power 
law at intermediate ( 14.5 < logA^Hi 5 , 16, iPetitiean et alj 
IT99I iHu et alJ IT995I: iKjmetal. 2002) and at flie highest 
column densities (log A^Hi 19, Storrie-Lomb ardi & Wolfd 
120001: iProchaska et al . 2005; O'Meara et al. 2007). A careful 
treatment of these systems is necessary, since even the inter- 
mediate column densities have a strong impact on th e total 
LyC absorption iMadaulll995l ; lHaardt & Madaij|[T996h . How- 
ever, due to the scarcity of 14.5 < logA^ni ^ 17 systems, the 
shape of the CD DF in this imp ortant range is presently not 
well constrained dKim et alj|200 2). 

We took a novel approach to constrain the high-z CDDF at 
14.5 < logA^Hi < 19 by matching the mean free path (MFP) to 
Lyman limit photons calculated from the CDD F to our recent 
measurements from SDSS at 3 . 6 < z < 4.2 (Prochaskaet al] 
I2OO9I see also [Prochaska et al.ll2010l) . The effective optical 
depth to L yman limit photons emitted at Zem and observed at 
Zobs is (e.g. lParesce et al.iri980l) 



feff ,LL (Zobs ,Zem)= / f (NuhZ) 



1 -e 



cWHidz,(l) 



with the Lyman limit photoionization cross section cJll = 
6.33 X 10^'^cm^ and the frequency distribution of absorbers 



in redshift and column density / (A^hi , z) = 



Consider- 



ing the different power-law redshift distributions of different 
absorber populations as outlined above, we approximated the 
CDDF as piecewise power laws that do not change over the 
considered redshift range, yielding 



feff,LL (Zobs , Zem) = ^Cj 



1 -e 



A'Hl.min.j 



(i+zym^ 



HI 



dA?Hidz,(2) 
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with different normalizations C/ and power law ex- 
ponents (YjiPj) in different column density ranges 
[Mii,min.>,-^Hi,maxj] • The normalization constants Cj 
are the products of the line density normalizations Aj 
(d«/dz = Aj (1 and the CDDF normalizations to yield 

an integral of unity in the respective column density range 



■'"HLmaxj' "HLmin.y 



(3) 



For the SLLSs we assumed jSsLLS = 1-4 ("O'Mea ra et alj 

12007). whereas for DLAs we adopted jSuLA = 2 

(iProchaska et alJ l2005l) . We fixed the contributions of 
SLLSs and DLAs to Teg^L with our explicit line density 
evolutions. These absorbers are highly optically thick to LyC 
photons, so their incidence rather than their column density 
distribution determines their share to TeffxL- 

By definition the MFP corresponds to the proper distance 
where Teff ll = 1 for Lyman limit photons emitted at Zem- 
In order to constrain the shape of the CDDF of Lyman 
limit systems and the Lya forest, we considered a contigu- 
ous triple power law at 12 < logA^Hi < 19 that results in a 
quasi-continuous CDDF over the full column density range. 
Requiring the 12 < logA^Hi < 19 forest to result in our as- 
sumed power-law redshift evolution of the effective Lya op- 
tical depth, we fixed y = 2.04 (7 = 0.20) for logA^ni < 19 at 
z > 1 .5 (z < 1 .5). We then varied the triple power law CDDF, 
each time simulating 1000 MC sightlines at < z < 4.6 in 
order to determine the normalization constants for the line 
densities Aj followed by computing the resulting total Ly- 
man limit effective optical depth (eq. |2]i and comparing the 
corresponding MFP a t 3.6 < z < 4.2 to our measurements 
(IProchaska et al.ll2009h . 

In order to find the most plausible values for the slopes and 
breaks in the CDDF we considered additional observational 
constraints. The CDDF is best determined in the z ~ 3 Lya 
forest and w e adopted i3i = 1.5 for 12 < logA^ni < 14.5 at 
z > 1.5 (e.g. lHu etalJIT995h . At logA^HL ,max,i — 14.5 we im- 
posed the first break in the CDDF to account for the d eficit 
of absorbers at logA^Hi > 14.5 (e.g. lPetitiean et al.lll993 f). Ini- 
tially, we tried a single )3 = 1 .5 power law that strongly under- 
predicted the MFP, but remarkably extrapolates into the SLLS 
and DLA range where the CDDF was set independently. This 
probably reflects the fact that the j3 = 1.5 power law approx- 
imation relies on both ends of the CDDF, which are by far 
the best constrained. We then varied the second break column 
density A^Hi.max,2 and the slope ^2 between the two breaks, re- 
quiring the slope jSj, at log.;VHi.max,2 < logA^Hi < 19 to meet 
the extrapolated j8 = 1.5 power law at logA^ni = 19, thus 
yielding a quasi-continuous CDDF, which we used at z > 1 .5. 

Our calculations confirmed previous results that the MFP, 
and thus the mean LyC absorption at high redshift, is very 
sensitive to the shage_of the CDDF at intermediate column 
densities (e.g. Madau 1995). In particular, we could rule out 
many parameter combinations (logA^Hi.max,2, j32) by requiring 
the calculated MFP to be consistent with both the normal- 
ization and the redshift evolution of the measured MFP at 
z > 3.6. In Fig.|5]we show our best match to the actual obser- 
vations, obtained for (logA^Hi,max.2i = (17.5,1.8), which 
imply a remarkably flat jSs ~ 0.9. The modeled MFP agrees 
extremely well with the observed values and can be accu- 
rately described by a power law at z > 3, yielding Amfp ~ 

50.14 [(1 +z) /4.6]^'^'^^ proper Mpc for a flat cosmology with 
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Figure 5. Compaiison of the mean free path /l„,fp as a function of red- 
shift z resulting from our adopted redshift evolution and column density 
distribution of Lya absorbers (open circles), direct measurements from 
Prochaska et al. (2009, filled circles) and the previous theoretical estimate 
by M adau et al. 1 1999) . All values are reported for a flat cosmological model 
with (£2m,iiA) = (0.3,0.7) and//o =V2kms"'Mpc"'. The solid line shows 
a power law fit to our empirical estimates /l„,fp = 50.14[(l+z) /4.6] 
Mpc, valid for z > 3. The dotted line shows the mean free path implied by 
a single /3 = 1.5 power law at logWni < 19 instead of the adopted broken 
power law. 



(£2m,ilA) = (0.3,0.7) and Hq = 72kms-' Mpc"'. In con- 
trast, by adopting a featureless P — 1.5 power law at log A^hi < 
19 together with the the slightly different distributions for the 
higher column density systems, the MFP is smaller by a fac- 
tor ~ 2.3 and is strongly inconsistent with the MFP measure- 
ments. The very good agreement between this underestimate 
and the MFP adopted by Madau et al. (1999) is not too sur- 
prising, as they assumed a single j3 = 1.5 and a single ab- 
sorber population evolving with redshift at 7 = 2, very similar 
to the 7 = 2.04 we adopted for logA^Hi < 19. We empha- 
size that at least two inflections in the CDDF are required 
at logA^Hi < 19 in order to yield a quasi-continuous CDDF 
that is consistent with our direct MFP measurements (see also 
Prochaska et al. 2010). 

Figure |6] shows the corresponding model CDDFs at z = 4 
(covered by our MFP measurements) and z = 2 (extrapolated 
from hi gher re dshifts using the redshift evolution laws from 
Section [3X2] |. The CDDF at logA^Hi > 19 is remarkably 
smooth, given that independent and uncertain redshift evolu- 
tion laws set the CDDF normalization there. The requirement 
for the CDDF to match the j3 = 1.5 power law extrapolation 
from the low column density forest at logA^Hi — 19 yields a 
continuous CDDF, both at z = 4 and at z = 2, as intended. 

As a consistency check we used our chosen distribution 
parameters to predict the incidence of Lyman limit systems 
(LLSs; logA^Hi > 17.2). Figure [T] compares the mock dif- 
ferential number densities of LLSs to observation s based on 
line counting ( Stengler-Larrea et al. 1995; Peroux et alj|2003l : 
iProchaskaet al.ii2010i: iSongaila & CowiQi201 0). Given the 
large statistical and systematic uncertainties in the observa- 
tions, the agreement is remarkable, even at 1.5 < z < 3.6, 
where we rely on the CDDF extrapolation from higher red- 
shifts. While our MFP measurements tightly constrain the in- 
cidence of LLSs at z > 3.6, the extrapolated CDDF might un- 
derestimate the incidence of LLSs if the CDDF straightens at 
lower z. By the same token, if LLSs evolve as strongly as in- 
dicated by Prochaska et al. (2010), we might have underesti- 
mated the MFP at z < 3.6. Our predictio n for the evolution of 
LLSs is most consistent with the fit bv iStengler-Larrea et alj 
( 119951) . who sampled z 3 based on earUer studies. Better 
models and predictions hinge on measurements of the MFP 
and the incidence of LLSs at z — 2-3. 
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Monte Carlo simulation parameters 



Figure 6. Modeled column density distribution functions (CDDFs) 
inj (d^HiAz) as a function of H I column density A'hi in a range Az = 0.1 
around z = \,Z = 1 and z = 4. For clarity the CDDFs have been scaled by the 
indicated factors. The symbols represent the binned CDDFs recovered from 
4000 MC sightHnes. The thick dashed lines show fits to the column density 
distribution at = 4 with the different slopes adopted in different column 
density regions (vertical dotted lines) designed to yield the measured MFP 
and its redshift evolution. The thick dotted line shows the CDDF at ^ = 4 
adopted for the forest (/5 = 1 .5) extrapolated to high column densities. Note 
that the redshift distributions and CDDFs in the SLLS and DLA range are set 
independently, whereas continuity is required for logA'm < 19. 
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'' The redshift evolution is parameterized by the effective optical depth 
Tcff .a = B ( 1 + z) ''^ ' or the line density dn/dz = A{1 + z)^ ■ 
^ The CDDF is a piecewise continuous power law dn/dNm « Wjjf ■ 
The b value distribution is dn/db = b^^ exp (— . 



umn density forest and slightly o verpredicts the frac tion of 
logA^Hi > 14.5 lines (their Fig. 5). iLehner et al] (I2007h found 
that the z < 0.4 CDDF steepens further at low column den- 
sities, whereas logA^Hi ^ 14.5 lines show a flatter slope 
j3 ^ 1.5. The low-redshift observations are inconsistent with 
our high-z model CDDF with its inferred low abundance of 
14.5 < logA^Hi ^ 17.5 absorbers. For simplicity, we there- 
fore assumed a featureless power law at z < 1.5 for the col- 
umn density range 12 < logA^Hi < 19, the slope of which 
was constrained by requiring a rough match to the CDDF 
at logA^Hi ~ 19 (set independently by the SLLS distribution 
from above), while yielding the observed number of LLSs at 
low redshifts jStengler-Larrea et alj|1995 j) and preserving the 
continuity in dn/dz for LLSs predicted from our extrapola- 
tion from higher redshifts (Fig.|7]l. A slope )3 = 1.55 matched 
these requirements. As an example, we show the modeled 
z = 1 CDDF in Fig. |6] A low er incidence of LLSs at z ~ 1 
as recently indicated by Songa ila & Cowid ( |2010 ) would not 
drastically change our predictions, because the total LyC ab- 
sorption at the He II edge primarily depends on the sparsely 
sampled redshift range z — 2-3. 

With our final set of input parameters (Table ^ we com- 
puted 4000 MC line lists over the relevant redshift range 
< z < 4.6. The number of sightlines is large enough to reach 
convergence in the incidence of optically thick H I absorbers 
even at low redshifts (Figs. |4]&|2]i, thus providing sufficient 
statistics for the highly stochastic UV LyC absorption. 



Figure 7. Differential number density dn / dz for Lyman Hmit systems 
(logA^Hi > 17.2) as a function of redshift z predicted from our MC 
simulations (open circles) comp ared to powe r law fits to actual data 
from Stengler-Lan'ea et al] j 19951) (dashed), Perou x et alj jlOOJi) (solid) and 
ISongaila & Co wie (2010) (thick dotted) in their quoted validity range. 
The filled sy mbols show a ctual mea surements I Stengler-Larrea et al. 1995; 
IProchaska'et al. 2010; Songaila & Cowie 2010). Note that the fit from 
[Songaila & Cowie (2010) includes previous data from Stengler-Larrea et al. 
Q995) at low z and Peroux et al. (2003) at high z- The vertical dotted lines 
mark the redshift range of our MFP measurements which constrain the num- 
ber of Lyman limit systems in the MC simulations. 



At low redshifts the CDDF is considerably less constrained, 
as the declining line density requ ires large samples of si ght- 
lines to be observed from space. iJanknecht et al.l (120061) de- 
termined a single power law for the CDDF at 0.5 < z < 1.9 
with /3 = 1.6, but their fit is dominated by the low col- 



3.2. Mock quasar photometry 

We used another IVlonte Carlo routine to generate mock 
quasar catalogs, i.e. distributions in emission redshift and ob- 
served magnitude, from the observed luminosity function of 
quasars. Due to the strong attenuation by the IGM, only 
quasars that are intrinsically bright in the continuum red- 
ward of H I Lya can be detected with current UV instru- 
ments. Thu s, we adopted the SDSS DR3 luminosity function 
( [Richards et al. 2006) that is well determined at bright magni- 
tudes. We integrated their Zem > 2.4 pure luminosity evolution 
model of the differential luminosity function in the observed 



/ band at redshift two 



combined with the co- 



moving volume in their adopted cosmological model to deter- 
mine the all-sky surface counts of quasars in a given range of 
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redshift and absolute magnitude 

(4) 

We chose to convert from absolute magnitude M;"'" 
to OT1450, the observed AB magnitude at 1450A in the 
quasar rest frame, vi a the relation M;"'"^ = M1450 — 1.486 
dRichards et al.ll2006l) . yielding 

mi450 =M;™=V51og(^^^ - 2.5 log (l+zem)+ 26.486 

(5) 

with the luminosity distance dL{zem)- By varying the in- 
tegration limits of Equation |4] we obtained a parameteriza- 
tion for dC4ji/dzem and dCiji/dmu^Q which we used to sim- 
ulate ^ 200000 pairs (zemj'WMSo) at 2.6 < Zem < 4.6 and 
15 < ;«i45o < 19. This large mock sample ensured an ac- 
curate sampling of the rare UV-br ight population of quasars 
transparent at the He II edge (see ! 14.1l below). For compari- 
son, Equation|4]predicts just ^ 11000 quasars on the full sky 
over the same range in redshift and magnitude. 

For each simulated quasar we assumed a unique spec- 
tral energy distribution mod eled as a power l aw /v °^ V^" 
with a break at H I Lya dTelfer et al.l 120021) . normalized 
to yield the modeled OT1450. Redward of the break we 
assumed a Gaussian distribution of spectral slopes with 
((oCcont) jCCofcont)) — (0.5,0.3) whereas blueward of the 
break we assumed ((auv) , c (ocuv)) — (1.6,0.6) consistent 
with the large variation in far-UV spectral slopes found by 
iTelfer et aTl (l2002:fl To the quasar continua we added the 
major quasar emission lines in the spectral range of in- 
terest (Lyj3, Lya, N V, Si IVh-O IV], C IV, C III], Mg II). 
The emission lines were modeled as Gaussian profiles, the 
strengths and widths of w hich were chosen consistent with 
IVanden Berk et al] (1200 Ih . with small variations from quasar 
to quasar. 

Lastly we blanketed each spectrum blueward of H I Lya by 
H I absorption in the IGM. For a given model quasar at a red- 
shift Zem we randomly drew one of our 4000 MC sightlines 
and computed the H I Lyman series and continuum absorp- 
tion at < z < Zem (§ 13.1.1b . yielding a final mock quasar 
spectrum at 912A< A <12000A. Blueward of He II Lya we 
assumed a He II Gunn-Peterson trough, resulting in zero flux 
(a reasonable assumption since the GALEX FUV band covers 
the He II break at Zem > 3.44). We th en obtained rnock S PSS 
ugriz photometry (asinh magnitudes. iLupton et al.l[T999h and 
mock GALEX FUV & NUV ph otometry (AB magnit udes) 
using the published filter curves jMorrissev et alj|2005l) . As 
Galactic extinction becomes important in the UV, we also 
computed the magnitudes after reddening each spectrum by 
the Galactic extinction curve (Carde niet ai]ll989i) . adopt- 
ing Rv =3.1 and a lognormal distribution in E{B — V) that 
closely r esembles the color exc ess distribution towards SDSS 
quasars dSchneider et al.l l2007t ). At the high Galactic lati- 
tudes considered here, the average extinction is ~ 0.26 mag 
and ~ 0.22 mag in the FUV and NUV, respectively. 

4. RESULTS 

' Note that lTelfer et al.l j2002l) quote the standai'd error of their mean spec- 
tral index instead of the (larger) standard deviation of the distribution of spec- 
tral indices (their Fig. 11). 
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Figure 8. Predicted cumulative probability distributions of the H 1 Lyman 
continuum optical depth TLyC at He 11 Lya in the rest frame of a source at 
redshift Zcm- 

4. 1 . The expected number of UV-bright quasars at z> 2.1 

For each of our ~ 200000 model quasars we calculated the 
total H I Lyman continuum optical depth TLyc at He II Lya 
in the quasar rest frame, i.e. the accumulated H I attenuation 
by [0.33 (1 +Zem) — 1] < z < Zem absorbers. This quantity 
characterizes the transparency of a sightline to the onset of 
the He II absorption, irrespective of lower-redshift LLSs that 
might truncate the spectrum in the He II forest region. We also 
computed the AB magnitude of the quasar at He II Lya 

m304 = mi45o + 0.191acoiit+ l-506auv + 1.086TLyC, (6) 

which depends on the input quasar magnitude mi45o, the spec- 
tral slopes of the continuum blueward (auv) and redward 
(oicont) of H I Lya and TLyc- 

In Fig. [H] we plot the cumulative distribution function of 
^LyC from our 4000 MC sightlines for different quasar emis- 
sion redshifts. Our calculations indicate a very low prob- 
ability to encounter a sightline that is not highly attenu- 
ated at the He II edge, consistent with previous estimates 
(Picard & Jakobsen 1993; Jakobsen 1998). The accumulated 
continuum optical depth strongly increases with emission red- 
shift. While ~ 9% of all quasars at Zem = 3 should be 'trans- 
parent' (TLyC < 1). this fraction drops to ^ 1% at Zem = 4. 

The predicted number of high-z quasars detectable in the 
far UV primarily depends on the increasing opacity and the 
declining quasar space density at Zem ^ 3. Figure |9] shows 
the predicted cumulative all-sky number counts of ;«i45o < 19 
quasars in the GALEX FUV & NUV bands compared to their 
predicted m304 for various redshift ranges. These estimates 
have not been corrected for Galactic extinction, in particu- 
lar close to the Galactic plane. For w E{B — V) > 1 com- 
monly encountered at Galactic latitudes \b\ < 20°, the FUV 
extinction is > 8 mag, so that ^25% of the sky are effec- 
tively blocked for He II s tudies even if quasars are found in 
this 'Zone of Avoidance' (lHubblelll934l) . 

The SDSS luminosity function predicts ^ 9200 mi4so < 19 
quasars on the entire sky at 2.7 < Zem < 4.5 (eq. H]). More 
than 200 of these should have OT304 < 21, well within the ca- 
pabilities of HST. However, at 4 < Zem < 4.5 there should be 
just ^ 600 mi45o < 19 quasars on the whole sky, the sightlines 
of which encounter larger LyC attenuation, yielding just ^ 1 
quasar at m304 < 21. At these high redshifts cosmic variance 
has a strong impact on the real number counts. The same is 
true for the least-attenuated UV-brightest quasars that are lo- 
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Figure 9. Predicted cumulative all-sky number counts of ^11450 < 19 quasars 
as a function of limiting UV magnitude for various redshift ranges. The up- 
per (lower) panel compares the predicted source counts in the GALEX NUV 
(GALEX FUV) band (solid lines) to the ones infeiTed at He II Lyce in the 
source rest frame (shown dashed in both panels), respectively. The gray 
shaded regions indicate Iff errors in the source counts due to cosmic vari- 
ance. Note that these predictions include the attenuation by the IGM, but 
they have not been connected for Galactic extinction. 

cated at the lowest redshifts. In order to obtain accurate results 
both at the highest redshifts and the brightest UV magnitudes 
we had to simulate the large set of ^ 200000 quasars, corre- 
sponding to 20x the predicted all-sky number counts. 

The GALEX bands trace the small UV-transparent quasar 
population very well, but differently at different redshifts. At 
Zem < 3.5 the GALEX NUV band is not a good indicator for 
the flux at He II Lya because of the high probability to en- 
counter a Lyman limit break in the large wavelength range 
between the NUV band and the onset of He II absorption. As 
the FUV band is closer to the He II edge it is a more sensitive 
indicator of flux at He II Lya. At Zem > 3.44 the FUV band 
samples the He II edge and the presumed He II Gunn-Peterson 
trough. The He II Lya absorption progressively attenuates the 
FUV flux and likely causes FUV dropouts at Zem > 4. Only at 
Zem > 4 will the GALEX NUV flux indicate a hkely transpar- 
ent sightline. 

Figure [TO] further illustrates the importance of detected 
FUV flux to select promising sightlines for He II absorp- 
tion. We show the normalized H I Lyman series and Ly- 
man continuum transmission spectra of four representative 
mock sightlines from the onset of the He II Gunn-Peterson 
trough to H I Lya at the emission redshift Zem = 3.4. In 
the sightlines shown in the upper three panels the indicated 
optically thick H I absorbers truncate the spectra at the Ly- 
man limit, causing dropouts in the overplotted filter bands. 
Obviously, only Zem ^3.4 quasars detected in the GALEX 



FUV band will show a transparent sightline that has recovered 
from intervening LLS breaks. Even for the small subset of 
high-z quasars detected by GALEX, intervening low-redshift 
LLSs likely truncate the quasar flux between the two GALEX 
bands. Thus, in order to select transparent sightlines at a high 
success rate, FUV detections are required at least at Zem < 3.4 
where the FUV band still samples the quasar continuum red- 
ward of He II Lya. 

4.2. Far-UV color selection of probable He II sightlines 

Figure[TO]also illustrates that the GALEX UV color mpuv ^ 
mNuv can be used to select the most promising sightlines to 
discover He II absorption. Significantly red GALEX colors 
indicate low-z LLS breaks (3rd panel of Fig.fTOli between the 
FUV and the NUV band, whereas blue GALEX colors signal 
the recovery from a LLS break or the relatively unabsorbed 
hard quasar continuum. NUV-only detections indicate trans- 
parent sightlines only if the FUV band significantly covers the 
strong He II absorption, i.e. at very high redshift (zem ^ 4). 

We used our mock quasar photometry to determine the 
fraction of transparent sightlines (defined as the fraction of 
sightlines with TLyc < 1) as a function of redshift. In Fig- 
ure [TT] we plot the probability contours that a quasar de- 
tected by GALEX at a given color will show a total TLyc < 1 
along the line of sight at He II Lya in the quasar rest frame. 
The UV-optical colors wnuv — '«I450 and wpuv — wJuso just 
give modest hints whether the quasar will show flux at He II 
Lya. The NUV-optical color indicates a transparent sight- 
line just at the highest redshifts, but is otherwise quite in- 
sensitive due to the frequent low-z LLS breaks between the 
NUV band and the He II edge. Thus, at any redshift the least- 
absorbed quasars with the bluest wnuv ^ w^i450 colors are the 
most promising candidates to detect He II. The FUV-optical 
color OTpuv — '«i450 provides better constraints. Quasars at 
■Zem < 3.4 at a mpuv ~ 'Wi450 ^3.4 have a > 60% chance 
to show a low TLyc < 1 along the line of sight. At higher 
redshifts, the He II Gunn-Peterson trough reddens the FUV- 
optical color 

The UV color mp\jy — wnuv (right panel of Fig. fTTT i yields 
the most natural color-selection constraints. Any quasar de- 
tected in both GALEX bands at a rather blue UV color has 
a high chance to show flux at He II Lya. Unless the FUV 
fluxes get severely absorbed by He II at Zem ^ 4, the GALEX 
UV colors of transparent quasars should be similar to those 
of their unabsorbed spectral energy distributions, with the 
slightly bluer colors indicating the recovery from partial LLSs 
that result in a steeply rising flux towards the FUV due to 
the strong frequency dependence of the LyC cross-section. 
Quasars at mpuv — '«nuv ^ 2 are likely to show a LLS break 
at the blue end of the FUV band even if they are detected in 
the FUV. Very blue quasars below the lower 20% line in the 
right panel of Fig.[TT]are recovering from a Tll > 1 LLS break 
at z > 2 so that their flux rises steeply in both GALEX bands. 

With these estimates on the UV color range of quasars that 
show flux at He II Lya, we can estimate He II detection prob- 
abilities for the actual GALEX-detected Zem > 2.7 quasars. 
Figure [T2]compares the GALEX mpuv — 'Wnuv colors of our 
transparent (TLyc < 1) mock quasars to actual observations. 
We find that the UV colors of quasars having sightlines that 
are known to be transparent down to the onset of He II ab- 
sorption are similar to the simulated UV colors of transparent 
quasars. A posteriori, the blue UV colors of most known He II 
quasars indicate a high probability for transparency. Among 
the GALEX-detected quasars without further follow-up the 
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Wav elength [A] 

Figure 10. Four simulated normalized H I forest Lyman series and Lyman continuum absorption spectra of a source at Zem = 3.4 with overplotted filter bandpasses. 
The upper three panels show sightHnes with optically thick systems that result in photometric dropouts without recovery at the He II Lya edge (vertical line). In 
the sightline shown in the uppermost panel, an intervening SLLS prevents the spectrum from recovering from the first encountered LLS break. The sightline in 
the lowest panel does not have an intervening optically thick absorber, so that the background source is detectable in all photometric bands until the onset of the 
assumed He II Gunn-Peterson trough. 
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Figure 11. Predicted colors of quasars detectable with GALEX as a function of emission redshift Zem i^^fi' ^NUV — ™i45o; middle: mpuv ~ '"uso; right: 
"ipuv — "iNUv). The contours delineate the probability that a quasar at a given redshift Zsm detected by GALEX at a given color will be transparent at the He II 
304A edge (TLyC < 1). The dashed lines mark the colors of the mean adopted quasar spectral energy distribution ( i]3.2t ignoring intergalactic absorption. 



rare quasars at mpuv — winuv ^ 1 are the best candidates to 
search for flux at He II Lya. Our MC simulations indicate 
a probability of > 60% that a Zem ^3.5 quasar detected at 
< nipuv — wiNuv ^ 1 will show TLyc < 1- The slight offset 
between the simulated and the observed UV colors of trans- 
parent quasars could be due to a generally harder UV spec- 
tral energy distribution than assumed in the simulations (i.e. 
(auv) < 1-6) and/or a higher mean LyC absorption from a 
larger population of z < 2 LLSs. We suspect the latter is more 
likely given the poor existing constraints on the exact CDDF 
and the evolution of the MFP at z < 3.6 ( ^3.1.4l i. 

In contrast, quasars confirmed by HST follow-up to show 
zero flux at He II Lya are mostly redder in mpuv — »inuv 
than the UV-transparent population, consistent with our sim- 
ulations. Especially the high upper limits wpuv — '«nuv ^ 2 
correspond to significant detections in the NUV, but no for- 



mal detection in the FUV, signaling the cutoff by an optically 
thick LLS. The only opaque sightlines that remain insensitive 
to our UV color selection are the ones intercepted by a LLS 
just within the narrow range between the blue end of the FUV 
bandpass and He II Lya (e.g. PKS 1442+101 in Fig.©. This 
is reflected in our simulations by the broadening color con- 
tours towards lower redshifts. 

If the LLS is not optically thick then the flux can re- 
cover, but the quasar is of very limited scientific value be- 
ca use it is too fa int for follow-up at He II Lya (i.e. TLyc > 
1; ISvphers et al.l 2009a b have identified two such quasars). 
Moreo ver, th e two BALQSOs confirmed in the FUV by 
show 



Sypher s et alj show red GALEX UV colors, presumably 
due to their intrinsically redder spectral energy distributions 
and/or BAL troughs extending in the UV. While these quasars 
are interesting to study the BAL phenomenon, they are effec- 
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Figure 12. GALEX UV color mpuv — '«nliv "f Zem > 2.7 quasars in the GALEX GR4 source catalog. For quasars in the SDSS footprint without UV follow-up 
observations we only show those that do not have sufficiently blue neighboring objects out to 5" separation to avoid GALEX source confusion. Open circles show 
quasars detected in both GALEX bands, whereas aiTows indicate upper (lower) color limits for sources detected in the NUV (FUV) only. Filled circles indicate 
the UV colors of quasar sightlines spectroscopically confirmed to be transparent at He II Lya. Quasars showing low spectroscopic FIJV flux, either due to a 
Lyman limit break in the FUV or a recovery from a Lyman limit break in the NUV, are shown as open squares. Thick crosses mark quasars targeted with HST, 
but not confirmed to show flux at He II Lya, either because of optically thick Lyman limit breaks in the FUV spectral range (crossed squares) or no detectable 
FUV flux at all (thick x crosses). Quasars without significant FUV flux in follow-up HST images are indicated as well (thin x crosses). Some spectroscopically 
observed quasars exhibit BAL features (diamonds), and several previous follow-up observations seem to have been affected by GALEX source confusion (+ 
signs). Thick star symbols mark the 8 UV-bright quasars we selected for upcoming follow-up spectroscopy with HST/COS in Cycle 17. The thick lines show the 
probability contours that a quasar at a given UV color shows flux down to the onset of the He II absorption, based on our MC simulations. 



lively useless for investigating intergalactic He II absorption 
as one cannot distinguish IGM He II Gunn-Peterson troughs 
from potential BAL troughs. 

The UV color separates well between blue He Il-transparent 
quasars and red opaque ones, despite the low S/N near the 
GALEX detection limit. However, quasars just detected in 
one of the GALEX bands require further attention. FUV-only 
detected sightlines probably recover from a partial LLS break 
so that the low NUV flux is beyond the detection limit. Given 
that we just quote 1 a flux limits on NUV dropouts, the colors 
of the 6 very blue confirmed He II quasars could be similar 
to those of the other He II quasars. Likewise, the FUV flux 
of some transparent quasars detected just in the NUV should 
have been detected as well. Nevertheless, since significant 
NUV-only detections indicate opaque sightlines, such back- 
ground quasars should not be regarded as prime candidates 
for spectroscopic follow-up. Generally, we do not consider 
very low S/N< 2 detections in a single GALEX band to be 
real, whereas sources detected in both GALEX bands proba- 
bly are, as the GALEX pipeline performs the source detection 
independently before merging the catalogs (Morrissev et al. 
|2007). 

Moreover, quasars with nearby optical neighbors should be 
avoided, as they will be likely affected by GALEX source 
confusion due to the broad instrument PSF. Apart from 
GALEX-detected Zem > 2.7 quasars with HST follow-up. 
Fig. [12] shows only those sources which qualify for further 
investigation (non-BAL, no blue neighboring source in SDSS 
DR7 at separation < 5"). Given that the UV color is not well 



constrained at low S/N (a S/N> 3 in both bands corresponds 
to a(OTFuv — otnuv) < 0.51) we consider two subsets of these 
quasars as the most promising ones for further detections of 
He II: (i) those 52 which have been significantly detected in 
both bands at S/N> 3 and have mpuv — otnuv < 1, and (ii) the 
1 14 remaining quasars detected at S/N> 2 in the FUV band. 
These samples are presented in Tables [3] and IH respectively. 
We caution that several GALEX detections outside the SDSS 
DR7 footprint will correspond to confused GALEX sources 
20% if we adopt our estimate from SDSS). Likewise, a 
few quasars with red optical neighbors (flagged in Tables [3] 
and|4|) might be confused sources, since our neighbor classi- 
fication was based on the broadband SED shape from SDSS 
and GALEX. Seven quasars detected on DIS survey plates 
have been flagged as potentially affected by source confu- 
sion. The majority (90%) of the Zem > 2.78 quasars in Ta- 
bles [3] and Bl were previ ously suggested as candidate He II 
quasars bv iSvphers et alj ([2009 a) . All but one of the 13 addi- 
tional Zem > 2.78 quasars h ave GALEX counterp arts beyond 
the match radius adopted bv iSvphers et alj[2009al (3")- 

5. APPLICATION TO THE SLOAN DIGITAL SKY 
SURVEY 

5.1. Comparing UV-bright SDSS quasars to predictions 

With a homogeneous well-characterized large area quasar 
survey such as SDSS, we can compare our predicted num- 
ber counts of UV-bright quasars to actual observations af- 
ter accounting for several observational effects. First, the 
predicted all-sky number counts (Fig. [9|l were corrected for 
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Galacti c fo reground extinction as incorporated in our calcula- 
tions ( ^3.2l i. Then we accounted for the actual GALEX GR4 
sky coverage and depth. The exposure time varies signifi- 
cantly among tiles of a given GALEX imaging survey, ren- 
dering them inherently inhomogeneous. Therefore we used 
the GALEX instrument sensitivity (Morrissev et al. 2007) and 
the actual GR4 tile exposure times to calculate a 5(7 limiting 
magnitude for each tile. With an approximate area correction 
for th e overlapping circular GALEX tiles (e.g. lBudavari et al] 
l2009h we then calculated the GR4 sky coverage as a func- 
tion of limiting magnitude. We regard the S/N> 5 threshold 
as sufficient to avoid incompleteness in the GALEX source 
catalog, but we note that apart from general source counts 
dBianchi et all 120071) the repeatability and S/N stability of 
GALEX is not well established at its instrumental limit. 

Next, we accounted for the SDSS sky coverage. Consid- 
ering that GALEX GR4 covers almost the full sky at high 
Galactic latitude (I/?! > 20°), we avoided the cumbersome 
calculation of the ac tual overlapping are a of SDSS DR7 and 
GALEX GR4 (see iBudavari et al.ll2009! for an appHcation to 
DR6H-GR3), and adopted instead the SDSS Legacy spectro- 
scopic sky coverage of 8032 deg2. SEGUE fields were not 
taken into account, as they are mainly at low Galactic latitude 
and have a significantly smaller quasar targeting rate. Lastly, 
we corrected for the SDSS quasar selection efficiency to pre- 
dict the number of UV-bright SDSS quasars detectable with 
GALEX. As SDSS selects quasars primarily by color, we used 
the photometric SDSS selection function by Richards et al. 
(|2006) averaged at / < 19.1. The magnitude cut provides a 
homogeneous survey limit at Zem > 2.7 (SDSS selects z > 3 
quasar candidates at / < 20.2) and ensures that the selection 
function does not depend on magnitude. Moreover, it is well 
matched to the rest-frame magnitude limit we applied in our 
simulations (/ ^ mi45o < 19), as the / band covers the quasar 
continuum redward of Lya at the relevant redshifts. The dif- 
ferent bandpasses induce a slight redshift-dependent offset 
i — nii4so —0.1, but uncertainties in the K correction used 
to determine the quasar luminosity function are larger than 
this. 

From our sample of quasars we then selected only those 58 
which were targeted by SDSS, have / < 19.0 and have been 
detected by GALEX at S/N> 5 in the NUV band. If we ex- 
clude SDSS quasars with blue optical neighbors that could be 
cases of GALEX source confusion, this number reduces to 52. 
Figure [T3] compares the cumulative number counts of these 
NUV-detected quasars to the prediction base d on our IGM 
model and the SDSS selection function by iR ichards et al.l 
(|2P06, upper curve). Adopting the Richards et al. (20()6|) se- 
lection efficiency, the number of NUV-detectable quasars is a 
factor of ~ 2 larger than observed, even including potentially 
confused GALEX sources. The predicted number counts can 
only be lowered by increasing the LyC opacity in our IGM 
model or by decreasing the SDSS selection efficiency. The 
uncertainties on other model ingredients, such as the lumi- 
nosity function of bright quasars, the K correction, and the 
GALEXh-SDSS footprint corrections, are too small to create 
this discrepancy. Given that our IGM model fits the MFP mea- 
surements (Fig.|5]) and independently reproduces the observed 
redshift evolution of LLSs (Fig. |7]) we have focused here on 
systematic effects in the SDSS selection efficiency. 

5.2. A color-dependent SDSS color selection function 

Quasar selection by broadband colors is expected to be 
inefficient and highly model-dependent at z ~ 3, where 
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Figure 13. Compaiison of predicted and observed cumulative number counts 
of Ztm > 2.7 i < 19 SDSS DR7 quasars detectable with GALEX at S/N> 5 
in the NUV. The thin lines plot the predicted number counts using the SDSS 
color selection efficiency from Richards et al. (2006, dashed), and our color- 
dependent selection efficiency (solid), respectively. Shaded regions outline 
Poisson errors on the source counts. The thick solid lines show the observed 
cumulative number counts of SDSS DRV w/o potentially confused GALEX 
detections due to nearby blue optical neighbors. 

quasar colors are sim ilar to those of main-sequence stars (e.g 
IRichards etani2006h . We used the simulated SDSS photom- 
etry of our ~ 200000 Zem > 2.6 model quasars ( 0.21 ) to re- 
assess the SDSS quasar selection function. SDSS selects most 
quasar candidates as outliers from the stellar locus in multi- 
dimensional color space. Because this procedure depends 
on the photometric errors, we computed these by fitting the 

J ihotometric erro rs of observed Zem > 2.7 SDSS DR5 quasars 
Schneider et al.i 2007) as a function of magnitude. We as- 
sociated each mock SDSS magnitude m with the fitted mean 
photometric error a,,, without modifying the mock magnitude. 
Thus, we assume perfect SDSS photometry with a realistic 
mean error, which simplifies our further discussion, but will 
likely result in an overestimate of the selection efficiency due 
to photometric uncertainties near the SDSS survey limit, par- 
ticularly in the u band. Potential effects of asymmetric distri- 
bution functions of SDSS magnitudes and their errors at the 
survey limit are beyond the scope of this paper 

Gordon Richards kindly agreed to process our mock pho- 
tometry with th e final SDSS quasar target selection algorithm 
(IRichards et al.l [2002) that incorporates the imposed 10% 
follow-up targeting rate of quasars whose colors intersect 
the st ellar locus (the 'mid-z' inclusion box of Richards et al] 
l2002h . The result of that operation is a selection flag for each 
mock quasar indicating whether it would have been targeted 
under SDSS routine operations. We then computed average 
SDSS selection efficiency as the fraction of selected mock 
quasars in Azem = 0.05 bins. 

In Fig. [14] we compare our selection function to the one 
by Richards et al.. (2006) . Both selection functions are es- 
sentially unity at Zem ^3.6 where colors of quasars are suf- 
ficiently red because of IGM absorption to separate well from 
the stellar locus. In particular, high-redshift LLSs will re- 
sult in red u — g colors due to u band dropouts (Fig. [TO] i. At 
Zem ^ 3.5, however, there is a striking difference between the 
two selection functions. Our average selecti on efficiency at 
Zem — 3.2 is ~ 25% smaller than the one bv IRichards et al.l 
( 120061) . whereas at Zem — 2.7 it is a factor of ~ 4 higher, and 
is in better agreement with their upper limit based on the ex- 
pected smoothness of the luminosity function with redshift. 
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Figure 14. Average SDSS quasar color selection efficiency derived using our 
MC siinulations of SDSS photometry (solid line) coinpared to the previous 
estimate by Richards et al. (dashed line). The average selection efficiency 
is almost magnitude-independent for the bright < 19.1 quasars considered 
here. The open circle marks the lower limit on the z — 2.7 selection efficiency 
estimated by Richards et al. ( 2006) by requiring a smooth luminosity function 
at lower and higher redshifts. The thick dotted line shows the SDSS color 
selection efficiency of niNUV < 26 quasars. The vertical bar marks the start 
of the first redshift bin we r egarded as unbiased for the measurement of the 
MFP ffrochaska et aU2009t) . 



The main model ingredients affecting the colors, and thus 
the selection efficiency, are the quasar spectral energy distri- 
butions and the IGM, i.e. the LyC absorption. Apart from 
a larger spread in the power-law spectral index blueward of 
H I Lya, our parameters to model the i ntrinsic quasar spec - 
tra are very similar to the ones used bv iRichards et al.l (120061) . 
so these discrepancies must be due to different assumptions 
regarding the properties of the IGM that result in statistically 
different quasar colors. 

The selection efficiency of our model quasars critically de- 
pends on the u — g color. Figure[T5]compares the distribution 
of mock u— g quasar colors to observed / < 19.1 SDSS DR7 
quasars ('Schneider et al."2010), either selected based solely 
on the Richards et al. (2002) color selection criteria, or on 
their radio flux. SDSS targets radio-detected quasar candi- 
dates independently of color The color-selected quasars have 
significantly redder u — g colors than the radio-selected ones 
at all redshifts Zem > 2.7 allowing for such a comparison (see 
the inset in Fig. [TSl l. They are also redder than most of our 
simulated quasars at Zem ^ 3.5, whereas the radio-selected 
quasars fill the simulated range in u—g color We verified that 
most quasars with very red u — g colors outside the simulated 
range are BALQSOs that were not treated by our MC simula- 
tions (for the selection efficiency of BALQSOs see lAllen et alj 
l2OT0h . 

The characteristic shape of the simulat ed color distribution 
is due to the SDSS magnitude system (iLupton et al.lll999h 
that yields finite values even for zero or negative fluxes. At 
Zem > 3.4 the frequent LLSs result in u band dropouts with a 
finite u — 24.63 at zero flux as defined for SDSS. At higher 
redshifts the g band flux is progressively attenuated by the 
IGM, which results in an artificially blue u — g color if the u 
band flux is zero. The m — g colors of Zem > 3.4 quasars are not 
well determined as the u magnitude exceeds usually employed 
detection limits. In this regime, the u band flux is systemat- 
ically overestimated due to Eddington bias, so that the ob- 
served u—g colors are bluer than the simulated ones without 
Eddington bias. Considering that SDSS selects even fainter 
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Figure 15. Mock vs. observed u — g quasar color as a function of redshift. 
Small filled circles (large open circles) show color-selected (radio-selected) 
/ < 19.1 quasars from the SDSS DRV quasar catalog (Schneider et aj 
[30l0). All plotted quasai's have been selected by the final SDSS criteria 
(Richards et al. 2002). Quasars marginally detected in the u band (CT„ = 0.2 
corresponds to S/N~ 5) are plotted in light gray. The dashed lines show the 
10%, 50% and 90% percentiles of the h — g color distribution of the simulated 
quasars. The solid contours show our derived SDSS selection efficiency at a 
given redshift and u — g color At 3.2 < Zem 5; 3.6 blue (u — g <2) quasars 
are missed by the SDSS color criteria, in contrast to radio-loud quasars se- 
lected independent of color The lines in the inset compare the cumulative 
probability distributions of color-selected (thin) and radio-selected (thick) 
2.9 < Zcm < 3 SDSS quasars. The latter exhibit systematically bluer u — g 
color than the former. 



/ < 20.2 high-redshift candidates, systematic effects in their 
colors at the faint end of the survey may non-trivially alter 
the selection function. Such effects are best explored by pho- 
tometric analysis of simulated survey images (e.g. Hunt et al] 
2004; Glikman et al. 2010). 

The thick contours in Fig. [15] show the SDSS selection 
efficiency at a given quasar redshift as a function of the 
u—g color. At high redshifts (zem ^ 3.6) the large range 
in color with a high selection efficiency means that almost 
all simulated quasars are selected regardless of their u — g 
color However, at 3 < Zem ^3.5 the SDSS quasar targeting 
algorithm preferably selects red quasars and systematically 
misses blue ones. This color-dependent selection efficiency 
is in good agreement with the distribution of the observed 
color-selected SDSS quasars in Fig. [15] In particular, very 
few observed quasars have u — g < I at Zem > 3, and most 
Zem — 3.4 SDSS quasars have u — g > 2, leaving a promi- 
nent 'hole' in color space compared to our predictions. On 
the other hand, the radio-selected SDSS quasars still reside 
in the color range of low selection efficiency. Our simula- 
tions also recover inhomogeneities in the color selection of 
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Zem ^ 3 quasars. iRichards et alj (120021) define the 'mid-z' in- 
clusion box at 0.6 < u — g < 1.5 with a targeting rate limited 
to 10% due to overlap with the stellar locus. However, can- 
didates having u — g < 0.6 are always follo wed up (this is the 
UV-excess criterion of Ri chards et al.ll2002l) . Hence, there is a 
'cluster' of DR7 quasars at m — ^ ~ 0.6 selected by UV excess, 
whereas at 0.7 < u — g < 1 there are very few color-selected 
quasars. 

The color-dependent selection efficiency of SDSS is due to 
the difficulty to differentiate quasar colors from stellar col- 
ors. The blue quasars at 3 < Zem ^3.5 do not separate well 
from the stellar locus, hence they are preferentially missed 
by the SDSS color selection criteria. But ho w does this ex- 
plain t he difference in the selec tion functio ns? IRichards et al.l 
( 120061) used the IGM model bv lFanl(ll99 9') that results in sig- 
nificantly redder u — g colors and a high selection efficiency of 
Zem ^ 3 quasars and, therefore, in a higher predicted selection 
efficiency. An explicit color dependence of the selection effi- 
ciency complicates 'completeness' corrections of color-based 
quasar surveys, rendering the average selection functions of 
Fig.[T4]invalid. To illustrate this further, we plot in Fig.[l4lthe 
average selection function of simulated quasars with a mea- 
surable NUV flux (mNuv < 26 including attenuation by the 
IGM). GALEX NUV-detected quasars are unusually blue in 
u—g, and consequently largely missed by the SDSS color se- 
lection criteria. 

The inefficiency of SDSS to select high-redshift quasars 
with blue optic al colors (and l ikely NUV flux) naturally ex- 
plains why the IRichards et al.l ('2006) selection function sub- 
stantially overestimates the number counts of NUV-bright 
quasars in Fig.[T3] Applying instead our color-dependent se- 
lection function lowers the prediction by almost a factor of 
2. Unexpectedly, the predicted number counts are now in ex- 
cellent agreement with the observed ones. In total, we pre- 
dict ~ 50 SDSS quasars in the DR7 footprint that can be 
detected at S/N> 5 in the NUV, very close to the actual 52 
(58) with (without) flagging potential cases of source confu- 
sion. We predict slightly too many NUV-bright quasars at 
21 ;$ wiNUV ^ 22, which may be due to the assumptions re- 
garding the quasar UV spectral energy distribution or the LyC 
opacity (the MFP is extrapolated at z < 3.6, Fig.|5]i. 

5.3. The SDSS Lyman limit system bias revisited 

Both, the observed differences in u — g color of color- 
selected and radio-selected quasars and the good match of our 
strongly color-dependent selection function to observations, 
point to significant selection effects of SDSS, either regarding 
the quasars themselves, or the intergalactic absorption along 
their lines of sight. As all relevant spectral parameters of the 
model quasars and all IGM absorbers along their sightlines 
were saved in our MC simulations, we could explore both pos- 
sibilities by comparing the statistical properties of the full MC 
sample and the subsample fulfilling the SDSS color selection 
criteria. 

Indeed, we find that the median UV spectral index of SDSS- 
selected model quasars is larger at Zem < 3.5 (Fig. \T6[. The 
Gaussian distribution of spectral indices is well preserved, but 
the mean is shifted to higher values, yielding redder u—g col- 
ors. Due to the increasing LyC opacity with redshift (see be- 
low), this bias decreases with increasing redshift. At Zem — 2.7 
there is a sharp break in the UV spectral index distribution of 
quasars that would be selected by SDSS. This feature can be 
attributed to the inhomogeneities in the SDSS targeting rate 
in M — g color space (Fig. [T5]l. Blue u—g colors can be due 
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Figure 16. Median spectral index a blueward and redward of H I Lya as 
a function of redshift as adopted in the MC simulations (solid) and for the 
subset selected by the SDSS quasar targeting algorithm (dashed). 

to hard UV spectral energy distributions, and the different 
targeting rates may cause non-trivial changes in the overall 
appearance of SDSS quasar spectra (e.g. composite spectra) 
as a function of redshift. The continuum redward of Lya 
is not significantly biased considering our simple model as- 
sumptions. Nevertheless, the slight shift to a harder spectral 
index at Zem < 3.5 might indicate too stringent selection cri- 
teria in the other three SDSS colors. We conclude that SDSS 
preferentially selects 2.7 < Zem 3.5 quasars with red spectral 
energy distributions in the u and g band. 

Lyman series and continuum absorption should have an 
even stronger impact on the u — g color at these redshifts 
(Fig. [TOl ). Therefore, we computed the mean IGM Lyman 
series and continuum transmission at different emission red- 
shifts, both for the full sample of 4000 MC sightlines and 
for the subsample of sightlines towards quasars fulfilling 
the SDSS color selection criteria in a Azem = 0.02 window 
around the emission redshift of interest. The resulting average 
'Lyrn an valley transmission spectra' (e.g. M0ller & Jakobseq 
119901) are plotted in Fig. \T7\ The sample of model spec- 
tra is large enough to clearly show the sawtooth-like fea- 
tures of overlapping Lyman series absorption. After an initial 
drop due to beginning series absorption at z < Zem the trans- 
mission recovers, because high-order high-redshift absorp- 
tion overlaps with low-order low-redshift absorption that de- 
creases with decreasing redshift. Beyond Lye there is a quasi- 
continuous roll-off of the transmission until LyC absorption 
sets in. At Zem = 3.6 there is essentially no difference between 
the average transmission of the full MC sample and SDSS- 
selected sightlines (compare the solid and dashed curves). 
However, at lower redshifts, the average LyC transmission 
towards SDSS-selected model sightlines is much lower than 
for general sightlines from the MC sample. The on average 
stronger LyC absorption corresponds to an on average redder 
M — ^ color Quasars at these redshifts are still in the vicinity of 
the stellar locus and LLSs in their sightlines will significantly 
redden the u—g color, moving them away from the stellar 
locus so that they can be selected by broadband colors. On 
the other hand, quasars with little LyC absorption (e.g. with- 
out LLS) will have colors similar to main-sequence stars, and 
are preferentially missed by broadband color selection. Due 
to the rarity of LLSs, however, their excess towards SDSS 
quasars should not significantly bias the Lya forest effective 
optical depth. 

Figure [17] presents further evidence that SDSS preferen- 
tially selects sightlines with strong H I absorbers at Zem < 
3.6, an effect that plagued the interpretation of our previ- 
ous results on the MFP and the number density of LLSs. 
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Figure 17. Mean IGM Lyman series and continuum transmission for sources 
at Zcm e {3.0,3.2,3-4,3.6} from our 4000 MC simulations (solid) and for 
the subset selected by the SDSS quasar targeting algorithm (dashed). Each 
sawtooth curve is due to accumulating Lyman series absorption (first Lya, 
then Lyo;+Lyj3, etc.), whereas the exponential-like roll-off is associated with 
the LyC opacity. Overplotted are the SDSS u and g filter curves (dotted). 

In iProchaskaltaLl dlOOl we found a significant flattening 
of the MFP at z < 3.6 that coincides with an apparent over- 
abundance of LLSs (Prochaska et al. 2010). We were puzzled 
that observed Zem 3.5 SDSS quasars are significantly redder 
than their brethren at Zem ^3.6 and suspected that the SDSS 
color selection criteria had biased our measurements. The 
mock quasar spectra processed with the SDSS color selec- 
tion routines support our previous claims. The SDSS-selected 
model quasars turn redder mu— g towards lower redshifts, as 
only these red quasars are outliers from the stellar locus. At 
3 ;$ Zem ^3.5 SDSS is essentially a Lyman break survey, re- 
sulting in an overabundance of LLSs and an underestimated 
MFP if the analysis is based solely on quasars from SDSS. 
Redder UV spectral indices of the quasars can alleviate this 
LLS bias somewhat, but not entirely. 

6. CONCLUDING REMARKS 

We have correlated verified Zem > 2.7 quasars to GALEX 
photometry to reveal the rare high-redshift quasars whose far- 
UV fluxes are not extinguished by intervening Lyman limit 
systems, with the goal to establish a sample of UV-bright 
quasars that likely show intergalactic He II absorption. We 
have used the GALEX UV color mpuv — otnuv to cull the 
most promising targets for follow-up. Red UV colors indicate 
that the quasar flux is prematurely truncated redward of the 
He II edge, whereas the rare quasars with blue UV colors and 
significant FUV flux will likely show flux at He II Lya. 

We have performed extensive Monte Carlo simulations to 
estimate the UV color distribution of UV-bright quasars and 
their surface density on the sky. We predict that ^ 600 200) 
quasars with Zem > 2.7 and muso < 19 should be detectable 
in the NUV (FUV) at muv < 21 (Fig.|9l without considering 
sources near the Galactic plane). The number of UV-bright 
quasars strongly declines with redshift due to the declining 
quasar space density and the increasing H I Lyman continuum 
absorption experienced at the He II edge (Fig. [8]l. Neverthe- 
less, there are enough targets within reach of HST/COS to 
significantly constrain He II reionization by He II absorption 
spectra, provided that the quasars are known and have been 
imaged with GALEX for efficient pre-selection. 

Most confirmed Zem < 3.5 He II quasars have blue UV col- 
ors and our simulations indicate a ^ 60% He II detection 



rate of quasars at similar UV color (Fig. [T2] |. a ^-^50% in- 
crease over approache s that do not include color information 
dSvphers et al.ll2009allB ). We have identified 166 additional 
quasars as prime targets for UV follow-up spectroscopy with 
HST/COS to significantly extend the sample of He II sight- 
lines before the end of HST's mission (Tables |3] and |4]i. We 
have started a survey with HST/COS in Cycle 17 to obtain 
FUV follow-up spectra of 8 UV-bright quasars selected from 
the much smaller GALEX GR3 footprint (star symbols in 
Fig. [nil. 

We have reassessed the SDSS color selection efficiency by 
applying the SDSS quasar selection criteria to mock photom- 
etry of our Monte Carlo spectra. We find that SDSS pref- 
erentially misses UV-bright quasars due to their blue colors 
that make them indistinguishable from main-sequence stars 
(Figs.fT4land[T5]l. The observed m colors of color-selected 
SDSS quasars are significantly redder than those of radio- 
selected ones at 3 < Zem ^ 3.5, and agree well with our color- 
dependent SDSS selection function (Fig. [TsT l. These missing 
quasars lack strong Lyman continuum absorption due to Ly- 
man limit systems along their lines of sight that would redden 
the u— g color (Fig.fTTIl. 

The SDSS color bias has not bee n well studied previously. 
Figure 18 of lBernardi et all (120031) reveals that SDSS rai-ely 
selected blue quasars at 3.2 < Zem < 3.6, bu t the a uthors did 
not investigate this further. [Richards et al.l (120061) explored 
whether primarily radio-selected SDSS quasars have differ- 
ent color selection efficiencies, but due to low number statis- 
tics, they regarded the differences to be insignificant (their 
Fig. 10). For the first time we have been able to demonstrate 
the full effect and its consequences. 

Since the UV-brightest quasars are among the bluest in 
SDSS u — g at all epochs, we conclude that SDSS is ineffi- 
cient in finding further promising targets for detecting inter- 
galactic helium. Although about two dozen quasars in the 
SDSS database alread y have been confirmed to show He II 
dSvphers et al.l l2009allbh . we predict that the FUV-brightest 
quasars without strong Lyman continuum absorption are in- 
sensitive to standard color selection techniques. 

Due to the restrictive SDSS selection criteria the statistics of 
high-column density IGM absorbers measure d towards color- 
selected SD SS quasars will be biased high dProchaska et al.l 
120091 120101) . Our results also indi cate that the incidenc e of 
DLAs based on SDSS samples (e.g iProchaska et al.l2005 *) has 
been overestimated. At 3 < Zem ^3.5 the few radio-selected 
quasars are probably the only ones within SDSS that are truly 
unbiased in the statistics of high-column density absorbers. 

The Lyman limit system bias will also affect the frequency 
of metal absorption lines that primarily occur at moderate to 
high H I column densities. Because these absorbers trace the 
large-scale structure of the IGM, the Lyman limit system bias 
might also impact analyses of the clustering properties and 
the power spectrum of the Lya forest. Presumably all Zem ?J 3 
quasars that have been first detected in broadband color sur- 
veys are affected by such a bias to some sort, depending on 
the exact color selection criteria and the number and effective 
wavelengths of the employed filters. Because the abundance 
of optically thick absorbers is far less constrained in reality 
than in 'completeness' simulations like ours, determinations 
of the optical quasar luminosity function at 3 < Zem < 3.5 
should invoke a variety of IGM models which will result 
in large systematic uncertainties in the luminosity function. 
Our study indicates that results based on color-selected high- 
redshift quasar samples are not as easy to interpret as previ- 
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ously thought due to the intertwined demographics of strong 
IGM absorbers and their background candles. 
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Table 3 

Promising targets to search for He II (GALEX S/N> 3 and mpuv — "iNUV < !)■ 



Object 


a (J2000) 


5 (J2000) 








filter 


niFUV [AB] 


fJiNuv [AB] 


limit'' 


S/Npuv 


S/Nnuv 


neighbors'^ 


2QZJ141 1—0229 


14''11™24?63 


— 02°29'42' 


'6 


2.702 


19.42 




22.52 


21.76 





5.3 


9.1 





FIRST J1007+4003 


10i>07"H6?85 


+40°03'56' 


'3 


2.728 


19.11 




20.33 


20.65 





4.6 


5.3 





SDSS J2330+0001 


23''30™26?26 


+00°01'23' 


'9 


2.733 


20.04 




23.79 


24.79 





12.3 


5.5 


3 


CTS 0216 


02*' 16^23'' 10 

\JA^ \.\J ^-J . \.\J 


— 39°07'55' 


'0 


2.735 


17.90 


B 


19.77 


19.37 





5.2 


1 1.9 


2 


SDSS J0029+0019 


00''29"H2?91 


+00° 1 9'46' 


'6 


2.736 


18.66 




20.68 


20.70 





25.6 


28.7 


1 


SDSS 10142—0077 


nil>42™43^54 


— 00°27'54' 


'0 


2.737 


19.93 




22.85 


22.58 





4.1 


5.2 





SDSS J2358— 0032 


23i>58™07?79 


— 00°32'24' 


'5 


2.753 


19.14 




21.67 


21.31 





9.9 


16.7 





UM 682 


03i>lO"28?10 


— 19°09'43' 


'7 


2.756 


17.90 


y 


19.68 


19.65 





7.2 


11.6 


2 


PC 2204+0127 


22'"06"M6M9 

\J\J ^VJ. I. y 


+01°41'45' 


'7 


2.757 


19.07 




22.33 


21.43 





5.0 


12.0 


2 


PC 1640+4711 


16''41™25?86 


+47°05'45' 


'8 


2.770 


19.51 




23.20 


23.06 





4.0 


4.8 





SDSS J2324— 0005 


23''24"52''55 


— 00°05'15' 


'3 


2.779 


19.43 




22.31 


22.16 





31.1 


32.9 


3 


SDSS J1309— 0333 


13''09"'34'' 1 8 


— 03°33'18' 


'4 


2.781 


19.19 




22.55 


22.24 





5.6 


6.8 





SDSS J0809+3116 


08''09™12?68 


+31°16'02' 


'1 


2.796 


18.86 


f 


21.27 


20.97 





3.1 


6.7 





O 0207—398 


02''09™28?59 


— 39°39'39' 


'5 


2.805 


17.15 


V 


20.29 


21.07 





5.2 


4.2 


2 


LBQS 1216+1656 


12''19"20?40 


+ 16°39'29' 


'5 


2.818 


18.17 




20.80 


20.04 





3.5 


6.9 





SDSS J15 19+3609 


15''19™10?37 


+36°09'40' 


'5 


2.819 


18.70 




20.87 


20.59 





5.2 


9.5 





Q 2315-4230 


23*' 18™ 15? 10 


— 42°13'48' 


'0 


2.830 


20.00 


V 


21.94 


21.37 





11.6 


21.0 


2 


SDSS J1230— 0253 


12''30™53?16 


— 02°53'52' 


'0 


2.837 


18.92 




21.88 


23.30 





7.0 


3.3 





2QZ J2158-3037 


2l''58"29?66 


— 30°37'21' 


'6 


2.838 


20.35 


b, 


24.57 


24.61 


1 


3.9 


3.3 


3 


HS 1024+1849 


10''27"34? 13 


+ 18°34'27' 


'5 


2.840 


17.83 




19.97 


19.82 





5.3 


16.9 





SDSS J0141 + 1341 


0li>41™34?01 


+ 13°41'58' 


'9 


2.843 


19.52 


Y 


22.95 


22.76 





3.5 


5.8 





SDSS J2156+0037 


2li>56™04?18 


+00°37'42' 


'3 


2.844 


19.01 




21.90 


21.39 





9.2 


13.3 


1 


CSO 0806 


13''04"H1?99 


+29°53'48' 


'8 


2.850 


17.65 


Y 


20.46 


20.47 





15.9 


26.1 





SDSS J233 1+0036 


23''31"31?48 


+00°36'44' 


'4 


2.852 


19.53 




22.55 


21.96 





4.1 


5.3 





SBS 1602+576 


16''03™55?92 


+57°30'54' 


'4 


2.858 


17.33 




19.76 


19.08 





8.2 


16.1 





PMN J1404+0728 


14i>04™32?99 


+07°28'46' 


'9 


2.866 


18.87 




20.78 


21.29 





4.0 


3.7 





PC 0058+071 5 

1. V. \J\J U U 1 \J 1 


n|l>00™58MO 


+02°31'32' 


'0 


2.868 


18.91 


n 


21.32 


21.22 





6.4 


7.4 


2 


CTS 0347 


22'"05"36''26 


— 34°26'03' 


'9 


2.870 


18.70 


n 


20.89 


20.48 





5.4 


7.9 


2 


2QZ J0126— 3124 


0ll>26™00?17 


— 31°24'21' 


'5 


2.881 


20.43 


hi 


21.48 


22.84 





4.6 


3.0 


2 


SDSS J1626+3856 


16''26™12?99 


+38°56'27' 


'2 


2.882 


18.63 




21.23 


21.82 





3.2 


3.1 





SDSS J2342— 0042 


-jVJ. y\J 


— 00°42'32' 


'8 


2.885 


20.47 




24.06 


23.21 





4.6 


4.2 





SDSS J1410+4727 


14'"10™59?61 


+47°27'33' 


'3 


2.901 


19.38 




21.74 


21.31 





3.1 


5.2 





SDSS J1443+3546 


14i>43miis5g 


+35°46'46' 


'3 


2.941 


18.79 




20.96 


21.19 





4.3 


4.6 





RDS 477A 


10'"53™06?04 


+57°34'24' 


'6 


2.949 


20.47 




24.55 


24.51 





3.1 


3.6 


3 


SDSS J08 18+4908 


08''18™50?01 


+49°08'17' 


'0 


2.954 


18.52 




21.49 


21.46 





11.6 


22.6 





SDSS J1033+5406 


10''33™10?71 


+54°06'46' 


'8 


2.959 


19.27 


Y 


22.55 


23.13 





4.9 


4.4 





FIRST J1456-0218 


14li5gm4Qsgg 


-02°18'19' 


'4 


2.963 


19.53 


Y 


22.80 


22.07 





4.5 


7.4 





SDSS J0922+5321 


09''22"M7?83 


+53°21'46' 


'6 


3.000 


19.75 


Y 


22.76 


23.32 





4.1 


3.4 





SDSS J1657+3553 


16''57™51?68 


+35°53'18' 


'0 


3.005 


19.23 




23.35 


22.81 





5.2 


7.4 





SDSS J0905+3057 


QghQjniQgsgg 


+30°57'57' 


'3 


3.027 


17.37 


Y 


20.89 


21.62 





6.5 


5.5 





SDSS JllOl + 1053 


1 1''01™55?73 


+ 10°53'02' 


'3 


3.031 


18.97 




21.59 


21.94 





4.0 


4.0 





SDSS J1244+6201 


12'"44"56?98 


+62°01'43' 


'0 


3.057 


18.63 


r 


21.08 


21.33 





4.7 


7.8 





SDSS J1052+2543 


10'"52"54?49 


+25°43'03' 


'9 


3.062 


18.52 


r 


21.43 


20.96 





3.1 


3.9 





PC 2211+0119 


22i>14"'27?81 


+01°34'57' 


'3 


3.100 


19.10 


R 


22.26 


22.35 





5.4 


6.7 


2 


SDSS J1025+0452 


10'"25™09?63 


+04°52'46' 


'7 


3.244 


18.02 


r 


21.37 


21.72 





3.8 


3.1 





SDSS J0955+6842 


09''55"'54?30 


+68°42'01' 


'2 


3.269 


19.26 


r 


24.05 


24.06 





5.4 


6.1 


3 


SDSS J1220+4549 


12'"20"H7?06 


+45°49'41' 


'1 


3.293 


18.20 


r 


22.78 


22.83 





3.8 


5.7 





HS 0911+4809 


09''15™10?01 


+47°56'58' 


'7 


3.337 


17.84 


r 


20.53 


20.25 





5.3 


9.4 





SDSS J0054+0028 


00'"54'"01?48 


+00°28'47' 


'7 


3.413 


19.93 


r 


22.06 


21.79 





7.6 


8.6 


1 


CLASXS 449 


10'"34""58?01 


+57°50'46' 


'5 


3.430 


23.80 


R 


24.19 


23.63 





4.5 


6.1 





SDSS J1233+0941 


12'"33"'02?74 


+09°41'44' 


'2 


3.816 


20.36 


r 


23.66 


23.03 





3.7 


3.4 





CDFN 097 


12'"36'"12?93 


+62° 19' 29' 


'8 


3.938 


22.80 


R 


25.78 


24.97 





4.5 


5.3 






Note. — 39 of the 41 Zem > 2.78 quasars were previously suggested as candidate He II quasars bv lSvphers et alji2009af) . 
" SDSS r AB magnitude if filter is r, otherwise Vega magnitude in given filter. 
GALEX Hmit flag. 0: formal two-band detection, 1: Iff lower limit in mpuv, 2: Iff lower limit in hinuv 

Neighbor flag. 0: no SDSS source within r < 5" of the quasar, 1 : sufficiently red SDSS source within r < 5" of the quasar, 2: quasar not imaged in SDSS 
DR7, 3: potential source confusion (DIS detection) 
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Table 4 

Further quasars with potential FUV flux (GALEX S/Nfuv > 2)." 







O [JZ\)\)\)) 




»'opt 


nilei 


'"FUV [ABJ 


"iNUV [ABJ 


limit 


o/lNpuV 


s/jnnuv 


neighbors 


2QZ J0035-2837 


00^35"" 24?23 


-28°37'14;'7 


2.702 


20.73 


bj 


22.70 


22.27 





2.2 


2.6 


2 


2QZ J0258-2941 


02''58'"09?15 


-29°41'08;'7 


2.702 


20.12 


bj 


21.78 


22.43 





3.5 


2.8 


2 


SDSS J1039+3040 


10''39'"24?05 


+30M0'59;'5 


2.705 


19.99 


r 


21.96 


23.78 


2 


2.9 


0.7 





SDSS J1301-0038 


13''01'"47?88 


-00°38'17;'3 


2.705 


19.40 


r 


22.99 


24.47 


2 


2.4 








2QZ J2153-2719 


2l''53'"16?08 


-27°19'38;'6 


2.706 


20.04 


bj 


22.14 


21.59 





2.5 


3.8 


2 


2QZJ0203-3153 


02''03'"15?58 


-31°53'54;'9 


2.710 


20.65 


bj 


22.21 


20.95 


2 


2.5 


2.0 


2 


CTS 0538 


14''21'"01?60 


-23°07'32;'0 


2.710 


18.50 


R 


21.31 


21.35 





2.7 


3.0 


2 


SDSS J1407+2127 


14h07'"01?12 


+21°27'15;'9 


2.711 


18.35 


r 


21.86 


20.88 





2.2 


6.8 





SDSS J1325+0814 


13''25""17?85 


+08°14'08!'4 


2.715 


18.70 


r 


22.31 


22.69 





2.8 


2.5 





SDSS J()014-0112 


O0'"14'"43?69 


-01°12'06!'4 


2.717 


18.84 


r 


23.99 


21.76 





3.7 


12.3 





Q 0040-370 


00i>42'"43?93 


-36°47'4i;'5 


2.723 


17.85 


V 


21.30 


21.33 





2.8 


3.8 


2 


2QZ J0141-3209 


0l''41'"54?69 


-32°09'li:'6 


2.724 


20.03 


hj 


22.77 


22.20 


2 


2.3 


4.3 


2 


SDSS J 11 59+0222 


ll''59'"04?30 


+02°22'14!'l 


2.725 


19.11 


r 


22.40 


21.10 


1 


3.5 


10.7 





Q 1613+172 


16''15'"56?87 


+ 17°07'5l!'4 


2.729 


18.24 


r 


22.00 


22.78 





2.8 


2.5 





QSO J0059-3541 


00'"59™14?21 


-35°41'42!'l 


2.730 


18.04 


V 


22.24 


20.97 





6.3 


16.4 


2 


SDSS J 1026+2842 


10'"26™54?39 


+28°42'54!'5 


2.739 


19.68 


r 


22.62 


22.09 





2.2 


3.7 





HE 0151-4326 


0l'"53"27?20 


-43°11'38"0 


2.740 


17.19 


hj 


20.63 


19.30 





5.8 


15.5 


2 


2QZJ1 129+0134 


ll''29™57?65 


+01°34'16"0 


2.743 


19.74 


r 


22.35 


22.49 


1 


2.2 


2.7 


1 


2QZ J 1326+0042 


13''26'"22M1 


+00°42'37:'3 


2.743 


18.72 


r 


21.66 


21.65 





2.7 


3.4 





2QZJ0053-3140 


00''53"'30?68 


-31°40'18;'8 


2.751 


19.99 


bj 


23.67 


21.87 





2.1 


7.1 


2 


2QZ J0012-3131 


00''12™43?11 


-31°31'13;'6 


2.755 


19.84 


bj 


22.64 


22.35 


1 


2.2 


3.2 


2 


HELLAS 149 


20''44'"34?80 


-10°28'08;'0 


2.755 


17.79 


V 


21.34 


20.15 





3.7 


8.0 


2 


QSO J0056-4013 


00''56'"11?76 


-40°13'16;'2 


2.758 


18.10 


R 


21.58 


23.15 





2.7 


2.2 


2 


SDSS J1600+4033 


16''00™33?09 


+40°33'43!'9 


2.761 


19.20 


r 


22.28 


21.58 


1 


2.4 


6.1 





SDSS J0150-0825 


0l''50°'09?46 


-08°25'10;'8 


2.763 


18.96 


r 


23.64 


23.29 





2.7 


3.4 





SDSS J0809+0658 


08''09'"46?14 


+06°58'07;'9 


2.763 


20.04 


r 


24.17 


24.01 


1 


2.1 


2.4 





SDSS J1546+2315 


15''46'"59?33 


+23°15'47;'3 


2.777 


17.80 


r 


22.22 


22.65 


1 


2.3 


2.9 





2QZ J0034-3048 


00''34'"47?21 


-30°48'13!'5 


2.785 


19.97 


h.i 


22.59 


21.46 





2.7 


6.4 


2 


SDSS J1418+5858 


14'"18'"22?89 


+58°58'06!'4 


2.785 


17.78 


r 


21.74 


19.94 





2.3 


11.7 





LBQS 0041-2707 


00'"43™51?83 


-26°51'27"5 


2.786 


17.83 


V 


21.79 


22.01 





3.3 


3.0 


2 


2QZ J0044-3147 


()0'"44™05?04 


-31°47'04:'5 


2.789 


19.80 


bj 


22.66 


22.27 





2.7 


5.5 


2 


2QZ J2223-3131 


22'"23'"12?45 


-31°31'29"4 


2.792 


19.44 


bj 


22.07 


22.59 





2.8 


2.8 


2 


SDSS J0103+0026 


0l'"03"37?46 


+00°26'08!'2 


2.795 


20.35 


r 


23.54 


25.65 


2 


4.4 








2QZ J1428+0010 


14'"28"49:'85 


+00°10'40!'7 


2.807 


19.79 


r 


21.65 


23.60 


2 


2.2 








H 0853+1953 


08'"56™26?47 


+ 19°41'37:'7 


2.818 


18.74 


r 


23.31 


22.17 





3.8 


6.5 





SDSS J0225 +0048 


02l"25'"19?50 


+00°48'23!'6 


2.820 


20.54 


r 


24.82 


23.01 





2.4 


5.7 





SDSS J0030+()053 


00'"30™17?11 


+00°53'58"8 


2.831 


19.92 


r 


23.78 


24.34 





3.3 


2.4 





FIRST J0905+3555 


09''05™36?07 


+35°55'5i:'6 


2.839 


18.39 


r 


21.88 


21.60 





2.3 


3.2 





SDSS J1504-0008 


15''04"'25?53 


-00°08'03!'2 


2.840 


18.92 


r 


22.44 


23.94 


2 


3.9 


0.9 





2QZ J0024-3149 


00''24'"16?22 


-31°49'42;'9 


2.846 


20.24 


bj 


22.56 


22.36 





2.6 


3.0 


2 


UM658 


22''46'"52?66 


-22°03'09;'2 


2.852 


17.80 


V 


22.40 


21.80 





2.3 


2.5 


2 


SDSS J0034-0109 


00''34'"20?62 


-01°09'17;'3 


2.854 


20.24 


r 


23.73 


22.55 





3.6 


6.5 





SDSS J1309+2815 


13''09'"39?49 


+28°15'08;'0 


2.854 


18.99 


r 


21.71 


23.17 


2 


2.7 


0.1 





SDSS J1439+0421 


14''39'"48?06 


+04°21'12;'8 


2.857 


19.00 


r 


23.99 


24.16 





2.6 


2.2 





SDSS J1241+2719 


12''41'"40?98 


-27°19'27;'5 


2.862 


19.21 


r 


22.54 


23.86 


2 


2.0 


1.1 





SDSS J0039+1527 


00''39'"39?96 


+15°27'20;'3 


2.867 


19.14 


r 


23.04 


23.99 


2 


4.0 


1.9 





FIRST J1231+0102 


12''31™39?12 


+01°02'29!'3 


2.883 


18.33 


r 


22.67 


20.57 





3.9 


26.2 





SDSS Jl 154+4030 


ll''54"13:'87 


+40°30'00:'l 


2.893 


20.36 


r 


21.64 


23.27 


2 


2.4 


1.0 





SDSS J0130-0007 


0l''30'°43?41 


-00°07'35:'3 


2.894 


19.95 


r 


23.76 


23.57 





2.8 


2.2 





2QZ JOl 14-2719 


0l''14"ni9?16 


-27°19'12:'4 


2.896 


20.55 


bj 


22.64 


23.30 


2 


2.6 


2.1 


2 


SDSS J 1322+3955 


13''22'"59:^97 


+39°55'29:'9 


2.898 


18.35 


r 


22.19 


23.10 





3.4 


2.7 





SDSS J1427+0014 


14'"27'"09?81 


+00°14'50!'2 


2.908 


18.54 


r 


23.56 


23.32 





2.6 


2.8 





PKS 0246-231 


02'"48"'22?74 


-22°57'58;'2 


2.914 


20.00 


R 


22.15 


21.42 





2.3 


3.1 


2 


SDSS J1525+2207 


15l"25"'34?50 


+22°07'00:'7 


2.914 


19.12 


r 


21.92 


21.87 





2.6 


2.8 


1 


SDSS J1210+3509 


12l"10"40?36 


+35°09'11"3 


2.919 


19.87 


r 


22.59 


21.72 


1 


2.1 


5.1 





FIRST J0936+2927 


09'"36™43?51 


+29°27'13!'6 


2.926 


18.11 


r 


20.80 


20.59 





3.0 


4.7 





SDSS J0300-0749 


()3ho()m47^62 


-07°49'02"8 


2.939 


20.02 


r 


22.84 


21.73 





3.7 


4.3 





FIRST J1604+1645 


16'>04"'41?47 


+16°45'38!'3 


2.939 


16.68 


r 


21.01 


19.56 





4.7 


16.2 





FIRST Jl 159+4136 


ll''59'"47?10 


+41°36'59;'l 


2.944 


18.71 


r 


22.12 


21.93 





2.6 


4.2 





FIRST J1332+0805 


13''32'"18?55 


+08°05'48;'3 


2.947 


18.86 


r 


21.90 


23.73 


2 


2.6 





1 


SDSS J0905+4107 


09''05'"18?02 


+41°07'57;'6 


2.954 


19.70 


r 


22.56 


23.01 


1 


2.0 


1.8 


1 


QSO J1334+2801 


13''34'"36?63 


+28°01'4i;'5 


2.958 


19.17 


r 


22.17 


23.85 


2 


2.1 


0.3 





SDSS Jl 143+3017 


lli'43'"14?67 


+30°17'li;'8 


2.964 


18.89 


r 


21.66 


23.43 


2 


2.5 








SDSS J2039-0047 


20''39'"06?09 


-00°47'36;'6 


2.966 


19.46 


r 


23.40 


26.03 


2 


3.2 


0.5 





SDSS J1335+2230 


13''35'"03?67 


+22°30'52;'7 


2.972 


18.95 


r 


21.42 


21.43 





2.7 


3.2 





SDSS J1356+0556 


13''56'"20?83 


+05°56'19;'7 


2.973 


19.03 


r 


21.73 


21.77 





2.6 


2.4 
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Table 4 

Further quasars with potential FUV flux (GALEX S/Nfuv > 2)." 
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SnSS 71040+2446 


10''40"'03*fi2 


+74°46'S1"0 

\^ "TU J J, U 


3.012 


19.46 




22.25 


22.1 1 





2.3 


2.1 





*\n's*\ 108^8-1-4.019 


KJO JO JJ.UZ. 


-1-40° 1 iHW" 1 

^'T'U i Z U J . i 
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n^to jZ tZ. 
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+ 14°21'16"9 
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19 70 


Q 
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2 
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1 9 
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2 4 
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-l-lfi°90'l4"4 
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3 122 
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3 148 
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3 
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n^iz Jl J/,U 
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+35°35'20;'6 
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r 


22.38 
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SDSS J2313+1441 
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3.337 


19.75 


r 


23.17 
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19.10 


r 


22.09 
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1 
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R 


23.44 


22.31 
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22.7 


3 
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3.410 


19.10 


hj 


21.53 


20.93 





2.2 
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2 


SDSS J1339+0703 


13''39'"51?84 
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3.438 


20.28 


r 


22.26 


22.67 


2 


2.1 
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1 
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+52°13'18;'0 


3.605 


18.77 


r 


21.90 
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1 


2.4 


2.0 


1 


Q 1422+231 
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r 
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21.81 





2.0 
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1 
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r 


21.96 


23.91 


2 
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1 



Note. — 76 of the 87 Zaa > 2.78 quasars were previously suggested as candidate He II quasars bv lSvphers et alj (2009af) . 
^ The sources listed in Table[3]are not repeated here. 

SDSS r AB magnitude if filter is r, otherwise Vega magnitude in given filter. 
^ GALEX limit flag. 0: formal two-band detection, 1: Iff lower limit in mpuv, 2: IfJ lower limit in m^uv 

^ Neighbor flag. 0: no SDSS source within r < 5" of the quasar, 1: sufficiently red SDSS source within r < 5" of the quasar, 2: quasar not imaged in SDSS 
DR7, 3: potential source confusion (DIS detection) 



